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Abstract 



D ■ Using algebraic tools inspired by the study of nilpotent orbits in simple Lie 

(~| ' 

algebras, we obtain a large class of solutions describing interacting non-BPS black 

fvg ' holes in AA = 8 supergravity, which depend on 44 harmonic functions. For this pur- 

^ ■ pose, we consider a truncation E'g(g)/5'pc(8, IR) C £'§(8) /5pm* (16) of the non-linear 

Cf^ ' sigma model describing stationary solutions of the theory, which permits a reduc- 

^— ^ ' tion of algebraic computations to the multiplication of 27 by 27 matrices. The lift to 

^ I A/" = 8 supergravity is then carried out without loss of information by using a perti- 

nent representation of the moduli parametrizing E-ji>j\/ SU^{^) in terms of complex 
valued Hermitian matrices over the split octonions, which generalise the projective 
coordinates of exceptional special Kahler manifolds. We extract the electromag- 
^ \ netic charges, mass and angular momenta of the solutions, and exhibit the duality 

^ \ invariance of the black holes distance separations. We discuss in particular a new 

type of interaction which appears when interacting non-BPS black holes are not 
aligned. Finally we will explain the possible generalisations toward the description 
of the most general stationary black hole solutions of AA = 8 supergravity. 
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1 Introduction 

The property that the weakly coupled calculation of the Bekenstein-Hawking entropy 
of BPS black holes [1] seems to generalise to non-EPS extremal black holes [2, 3, 4], 
suggested an increasing interest in characterising the entire moduli space of such solutions 
in various supergravity theories. The understanding of non-BPS extremal black holes is 
an important step toward the study of black holes carrying a non-zero temperature. 
The classification of supersymmetric composite black hole solutions has permitted to 
understand the mismatch between the enumeration of spherically symmetric BPS black 
holes in A/" = 2 supergravity and the counting of BPS states within weakly coupled 
string theory [5, 6]. The associated wall crossing formula can also be understood from 
the moduli space of regular solutions [7, 8]. Because the weakly coupled computations 
for non-BPS black holes are not based on well established non-renormalisation theorems, 
it is important to extract the maximum of information from the classical solutions in the 
supergravity limit. These solutions have been studied extensively in the literature, see 
[9]-[34] for a non-exhaustive list of progress that have been achieved in the recent years. 

The approach we will follow relies on the well established result that the stationary 
solutions of 'irreducible' supergravity theories with abelian vector fields are described 
within a non-linear sigma model over a pseudo-Riemannian symmetric space G/K* cou- 
pled to Euclidean gravity in three dimensions [35]. Here we will only discuss A/" = 8 
supergravity, therefore we will fix notations accordingly and consider the coset space 
-E'8(8)/S'pm*(16); although note that the generalisation of the content of this paper to 
exceptional M = 2 supergravity theories will be completely straightforward in our no- 
tations. The scalar momentum P is defined as the component of the Maurer-Cartan 
form V~^dV in the coset component 128 = e8(8) 0so*(16) of the Lie algebra. Solu- 
tions describing spherically symmetric black holes are then determined by the associated 
Noether charge in the Lie algebra e8(8) of -Escs), and can therefore be classified in terms 
of -E8(8)-orbits [35, 36, 37]. In the extremal limit, the Noether charge is nilpotent and the 
spherically symmetric extremal black hole solutions are classified in terms of the class of 
nilpotent orbits of -^8(8) in ^8(8) which lie in the closure of the minimal semi-simple orbit 
of -Escs) in ^8(8) [37, 38]. Then P automatically lies in the same nilpotent orbit, and more 
precisely in its intersection with the coset component 128. 

Relying on these properties, it has been proposed that the most general stationary 
solutions which geometry is fibered over a fiat three-dimensional base are described by 
fields V valued in specific nilpotent subgroups [34], such that P is automatically nilpotent. 



Indeed, a vanishing three-dimensional Riemann tensor imphes the equation 

Tr P,P, = (1) 

and non-nilpotent solutions would then involve imaginary eigen values, which have been 
shown to lead to singularities in the spherically symmetric case [36]. Moreover, all known 
under-rotating extremal solutions (as opposed to over-rotating extremal solutions gener- 
alising the extremal Kerr black hole) satisfy to this criterion, as it has been shown for 
the BPS solutions in [39], and for both the almost BPS solutions introduced in [12] and 
the composite non-BPS ones in [34]. 

In this paper we will exploit these ideas to derive generalisations of these solvable 
systems of differential equations to compute and study a large class of solutions in A/" = 8 
supergravity. The generalisation to A/" = 8 is pertinent because this is very probably 
the simplest such supergravity theory at the quantum level, and on the other hand, 
it is complicated enough classically such that its solutions exhaust the possible classes 
of solutions one can have in supergravity coupled to abelian vector fields and scalars 
parametrizing a symmetric space. We will study the solvable systems of equations that 
are associated to the nilpotent orbits of C8(8) which can be realised within so (4, 4), and 
which are therefore the direct generalisations of the systems of equations studied in the 
STU model in [17, 14, 34] to the full M = 8 supergravity. However, as opposed to the BPS 
solutions, the non-BPS solutions do not in general sit in a given M = 2 truncation of the 
theory, and our solutions will include new harmonic functions. Although these functions 
will not dramatically change the physical properties of the solutions, they modify the 
systems of differential equations in a non-trivial manner. 

Within the maximal A/" = 2 truncation of AT = 8 supergravity, the relevant nilpo- 
tent subalgebras defining these solvable systems admit an SU{2) x^^ >S'L(3,H) automor- 
phism, and they decompose accordingly into a graded algebra; where 5'L(3, H) = SU*{Q) 
is the group linearly realised on the special coordinates of the special Kahler space 
S'0*(12)/f/(6) of the maximal M = 1 truncation of the theory. The components of 
these nilpotent subalgebras in the coset component (which are associated to the har- 
monic functions defining the solutions) are 
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whereas for A^ = 8 supergravity one has 

nA,_.BPs = (1 © 15)'^' © (2 ® 6)'^) © 15^'' © 1^^' 
n„o„-BPs = (2 © 15)'^> © (2 © 6)'^' © 2'^> 

Ubps = (2 © 6)(^' © (1 © 15)^^' © (1 © 15)<^' (3) 

such that SU{2) x^j SL{3,M) is also an automorphism of these nilpotent subalgebras. 
The integer superscripts indicate the grading preserved by these nilpotent algebras. It is 
such that the new generators of grade 2 are simply abelian in the generalised non-BPS 
nilpotent algebra. In this sense, the generalisation of the non-BPS system is the simplest. 
The almost EPS algebra is the next, with new generators of grade 2 that only commute to 
source the grade 5 scalar function. The BPS algebra is the most stringent generalisation, 
and the new generators of grade 1 then source half of the old ones. We will see that in the 
associated, say 'locally BPS system', these new generators are associated to axion fields, 
and would source strings instead of black holes. The more general nilpotent subalgebras, 
which will be discussed in the last section, always have the property that new generators 
appear at grade 1, and modify strongly the algebra. It is not yet clear if they can describe 
composite black hole solutions. 

We will make progress in the physical understanding of these solutions, which are also 
pertinent within the M = 2 truncations. In particular we will generalise the two non-BPS 
centres solution of [34] to an N centres solution. We will not assume the solution to be 
axisjTiimetric, although we will only be able to derive a closed form formula in the case 
in which the centres are all aligned. The missing explicit function will be defined as a 
convergent integral for which we will derive the asymptotic expansion in the near horizon 
and the asymptotic regions in Appendix A. Therefore we will not prove the regularity of 
the solution everywhere, but only in these different regions. We will nevertheless be able 
to prove the absence of closed time-like curve due to Dirac-Misner string singularities. 
These criteria are known to be enough for regularity by experience, but a complete proof 
would require a careful analysis of the solution everywhere. We can nevertheless rely 
on numerical simulations that have been carried out in [17, 34] to exhibit that these 
solutions are indeed generically regular everywhere. 

In the non-BPS composite case the local properties of the non-BPS black holes at 
the horizon are as expected. In the case in which they are not aligned, they produce 
nevertheless an additional angular momentum which is related to the interactions of the 
centres three by three. Studying the equations for the distance separations between the 
black holes, we will see that they are not linear in their inverse distance separations. 



such that the interactions between the black holes do not only depend on their distance 
separations pairwise, but also on the geometry of the triangles joining them in triplets. 

We will study the composite non-BPS solutions in much more details, and in par- 
ticular we will have an additional section devoted to the study of the duality invariance 
of their physical properties. This will permit to get some information on the general 
case, not associated to a specifically simple duality frame in which the solution can be 
written in closed form. In this section we will exhibit the general form of the momentum 
P G 128, which will permit to introduce a generalised fake superpotential for A/" = 8 
supergravity. The introduction of auxiliary tensors that can be determined in terms of 
the electromagnetic charges of the black holes and the asymptotic moduli will permit to 
show the duality invariance of the equations determining the distance separations of the 
centres. 

In the almost BPS system we will study the restricted case of one EPS centre and 
arbitrary many non-BPS ones. This solution will be complementary of the one derived in 
[14], where they considered one non-BPS centre and arbitrary many aligned BPS ones. 
We will see that even if this system does not allow for interactions between the non- 
BPS black holes in the absence of BPS black hole, the presence of the latter permits to 
consider charge configurations which do not locally commute for the various non-BPS 
black holes. The local properties of the BPS black hole in the near horizon region are 
somehow unexpected, because we will see that there is no enhancement of supersymmetry, 
and the BPS horizon is replaced by a surrounding under- rotating horizon,^ which area is 
larger than the one of the original BPS centre. In this case also we will find additional 
contributions to the angular momentum which appear when two non-BPS black holes 
are not aligned with the BPS centre. In general one should expect to have additional 
corrections when three BPS centres are not aligned with two non-BPS ones coming from 
the solution of the Laplace equation with five distinct point sources. 

In order to carry out this program we will not consider directly the non-linear sigma 
model defined on Es(^s^ / Spin* (16) . It turns out that the smallest linear representation 
of -Escs) is the adjoint which is 248-dimensional. Instead, we will consider the truncation 
to the Eq(^q^/SPc{8, R) coset space. This truncation contains all the information we need 
and the 27 dimensional fundamental representation of -^6(6) can be managed with the 
help of a computer. It will appear that the SL{3, R) linearly realised symmetry of 



^By under- rotating horizon we mean an ergo- free horizon with zero angular velocity, but which is not 
spherically symmetric such that it would lead to an ADM angular momentum for a single-centre black 
hole [40]. 



the solutions in the truncation is promoted to a hnearly reahsed SU(2) x SL{3, H) = 
SU(2) X SU*{6) symmetry in A/" = 8 supergravity. Relying on this property, we will be 
able to generalise the solutions to A/" = 8 supergravity without ambiguities. Doing so, 
we will define particularly convenient coordinates for the symmetric space E-jfj^/ SU^^) 
which generalise somehow the special coordinates in very special Kahler geometry. The 
basic idea is to realise Ej(^j)/SU^{8) as a Ej(^j)/{U{1) x -^6(2)) fibered over the quaternionic 
space SO{4,4:)/{SO{4:) x ^0(4)) 

E7(7)/(f/(l) X i?6(2)) ^ i?7(7)/5t/c(8) 

i (4) 

50(4,4)/(50(4) X 50(4)) 

The complex fibre £'7(7)/(?7(l) x i?6(2)) is a pseudo-Riemannian special Kahler space, and 
can be coordinatized in terms of Hermitian matrices over the split octonions, similarly 
as special coordinates for the exceptional special Kahler space Ej(^_25)/{U{1) x Eq(^^js))- 
However, the fibre-bundle structure does not preserve the complete £'7(7) isometry of the 
fibre. Therefore the associated special coordinates will not admit a linearly realised -Eece) 
symmetry preserving the equations of motion, as the special coordinates admit a linearly 
realised -E'6{-26) symmetry in the exceptional special Kahler geometry [41]. 

In this paper we will parametrize the scalar fields in terms of a complex exceptional 
Jordan algebra element t, i.e. a 3 by 3 Hermitian matrix over the split octonions, which 
we will decompose into a 3 by 3 Hermitian matrix over the quaternions tj and an anti- 
symmetric matrix over the quaternions t2 (equivalently a 3-vector of quaternions) 

t = ti + hz (5) 

and similarly for the electromagnetic charges Qq, Cl,'P,p^ and the electromagnetic fields. 
This notation will be rather convenient in order to render straightforward the truncation 
of our solutions to A/" = 2 supergravity. In A/" = 2 supergravity, the components tz linear 
in the split imaginary unit £ are simply set to zero. The generalisation to any M = 2 
supergravity coupled to vector multiplets which scalar fields parametrize a symmetric 
special Kahler space of cubic prepotential can be obtained by considering tj (and the 
other Jordan algebra elements) to lie in the appropriate Jordan algebra [41, 42]. In 
particular, one can obtain the solutions of the STU model by considering the Jordan 
algebra of diagonal 3 by 3 real matrices, such that the three diagonal components of t 
give the three STU moduli, and respectively for the charges. 

Similarly, the SL{Q)/ SO{Q) moduli of the truncated theory will be described in terms 

6 



of special coordinates associated to the fibration 

SL{6,R)/{U{l)x SL{3,C)) ^ SL{6,R)/SO{6) 

i (6) 

which fibre will be coordinatized by a complex element of the Jordan algebra of Hermitian 
3 by 3 matrices over the split complex. The solutions described in this paper will admit 
a linearly realised 5*^(3, R) symmetry, such that the two base space scalars are not 
sourced, and the solutions can be written in terms of symmetric 3 by 3 real matrices 
and a real 3 vector. The enhancement of this symmetry to SU{2) x SL{3, H) in A/" = 8 
supergravity is very constraining, and implies that the corresponding explicit solutions 
can be determined without ambiguities in terms of 3 by 3 Hermitian matrices over the 
quaternions, and a 3 vector of quaternions, such that the 16 S'0(4, 4)/(S'0(4) x 5*0(4)) 
base space scalars are not sourced. For the most general ergo-free extremal solutions 
of A/" = 8 supergravity discussed in the last section, the linearly realised symmetry is 
reduced to 5*0(4,4), which is a priori not big enough to determine without ambiguities 
the generalisation of a solution of the truncation to A/" = 8 supergravity. Nevertheless, 
combining the knowledge inherited from the most general solutions of the truncation with 
the explicit form of the first order linear system inN' = 8 supergravity, one might be able 
to determine the most general solutions without considering directly the £'8(8)/5'pm*(16) 
representative. 

2 Exceptional truncations 

In order to find solutions of A/" = 8 supergravity in four dimensions, and more gener- 
ally of eleven-dimensional supergravity, it is often extremely useful to consider pertinent 
consistent truncations of the theory. Because one has been mainly interested in BPS 
solutions preserving a certain amount of supersymmetry, one was used to consider super- 
symmetric truncations of the theory. However, we will be studying non-BPS solutions in 
this paper, and there is no particular reason to consider truncations which are themselves 
supersymmetric. On the contrary, it will turn out that the most appropriate consistent 
truncations of A^ = 8 supergravity which are pertinent in the discussion of composite 
black hole solutions do not define the bosonic sector of supersymmetric theories. Rather 
than maintaining a certain amount of supersymmetry, these theories preserve somehow 
"exceptionality" . 

7 
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A^ = 8 supergravity can be obtained as the compactification of eleven- dimensional 
supergravity on a seven-torus [43], or alternatively as the compactification of type IIB 
supergravity on a six-torus. The second theory admits a non-supersymmetric truncation 
which consists in setting to zero the scalar fields and the doublet of 2-form fields. The 
remaining theory consists in a self-dual 5-form field strength coupled to gravity. 

A further truncation can been carried out in eight dimensions, by considering gravity 
coupled to an SL{2)/SO(2) non-linear sigma model and a 3-form, which 4-form field 
strength transforms in the fundamental of SL{2) (by complex selfduality) . 

These theories have the common property that after dimensional reduction down to 
four dimensions, the corresponding stationary solutions admit a duality group of type 
En(n)- Similarly as the exceptional M = 2 supergravity theories, their respective duality 
groups in 3, 4 and 5 dimensions define the magic square 



(7) 



associated to the three exceptional Jordan algebras of Hermitian 3x3 matrices over 
the three spit composition algebras. The realisation of split real form of exceptional Lie 
algebra as quasiconformal groups associated to exceptional Jordan algebras defined over 
split composition algebras has been defined and studied in [44, 45]. 

The stationary solutions of the corresponding four-dimensional theories are described 
by three-dimensional Euclidean gravity coupled to the non-linear sigma model over the 
pseudo-Riemannian symmetric spaces 

i?6(6)/5pc(8,R) C i?7(7)/^f/c(4,4) C Es^s)/Spin:{16) (8) 

respectively (where the index c states for the quotient by the Z2 subgroup which leaves in- 
variant the real 27, the complexe 28 and the real 128 representations of 5*^(8, R), S'f/(4, 4) 
and S'pm* (16), respectively). 

The main advantage in considering these truncations is that the smallest linear rep- 
resentation of -£^8(8) is of dimension 248, whereas -Eece) admits a 27-dimensional linear 
representation. Moreover, one will be able to lift the solutions associated to the -Eece) 
model to the general solutions of A/" = 8 supergravity by promoting real functions to 
functions defined over the quaternions. 

For this purpose, let us first define the coordinates over Eq(^q)/Sp^{S, R). 



2.1 Conventions for £^6(6) 

In four dimensions, this theory consists of Einstein gravity coupled to 20 scalar fields 
parametrizing the symmetric space SL{6)/S0{Q) and 10 vector fields which field strengths 
transform altogether in the 3-form representation 20 of SL{Q). The stationary metric 
will take the standard form 



as well as the vector fields 

VsAabc = C,ABc{dt + cu) + wabc , 



(9) 



(10) 



where the tensors are antisymmetric with respect to the permutations of the three SL(6) 
indices ABC. 

We will decompose the 27 representation in terms of the four- dimensional duality 
group as 



e6(6) = 1'-'' © 20'-'^ © (gli © sky' © 20<^^ © 1 
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such that a general element of C6(6) acting on a vector (i?'^', S^ i)^ j^c)) admits the following 
matrix form in terms of two antisymmetric tensors Eabc and Fate, and the traceless tensor 

b 
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where Eabc = l^abcdefE'^'^-^ and respectively for Eabc- The elements of sp(8, R) are defined 
from the involution associated to the metric 



as X- 
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where /!"(, is antisymmetric. 

Using the standard convention that capital indices correspond to rigid SL{Q) (acting 
on the left) whereas small ones correspond to local 50(6) (acting on the right), one has 
the representative 



V = exp [('^''^Eabc + (tE] exp[f/H] 
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where one must note that Va"^ should be understood as the transverse w* acting on the 
right through the contraction with the index a. 

It will be convenient to make use of a non-triangular representative in SL{Q)/ SO{Q) 
as 
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which gives 
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where t* = a, + ze"^"^' are the STf/ moduh. The reason for this choice is that the matrix 
on the right corresponds to the GL{1) x SL{3)/SO{3) x SL{3)/SO{3) representative 
(with the two SL{3) factors respectively in upper and lower triangular form), which 
parametrizes the five-dimensional scalars and the dilaton (p = ^^ (pi that appears in the 
dimensional reduction from five to four dimensions, whereas the left matrix corresponds to 
the axion fields associated to the five dimensional electromagnetic fields, which transform 
in the 3 (S) 3 with respect to the five-dimensional duality group SL{3) x SL{3). 



Exceptional Jordan algebra 

This representation will permit to use the homomorphism SL{3) x 5*^(3) = SL{3, C*) 
where C* is the split complex composition algebra, with imaginary unite i satisfying 
i'^ = 1. In this way one can rewrite the axion fields as a (Jordan algebra) Hermitian 
matrix over C* 

' ai as 02 
a = at, a2 (X\ (18) 

\ ^2 a\ a-s J 

with 



ai = -{^o+e)+^{e-^o) 



a2 = -(C' + Co) + 2(C'-Co), as 



(19) 
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Using the property (which defines the homomorphism C* = R, © R) 

■1±£\2 1±£ 



(20) 



one computes that the components of the Maurer-Cartan form v^ ^dv^ associated to the 
axions take the form 
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e = ^(e^-6)-^(e^ + 6), C = ^(C3-C^) + ^(C3 + C^), v = \{v'-vz)-{{rf + ^,) 

(22) 
and p(^) acts on the resulting matrix such that on any Hermitian matrix 
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with 



p,(,) = 6-^(2^++-) , p2(0 = e^(-^++-) , p3(,)=e-^(^++2-) 



the three phases satisfying the triahty condition 
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The matrix 
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on the right of (21) also defines partially the GL{1) x SLi?)^ C*) coset representative. In 
analogy with the A/" = 2 case for which the GL{1) x SL{'i^ C) Maurer-Cartan form can 
be written as 

dvv^^ + v^^^dv^ = v^^^d{v^v)v~^ (27) 



we will define the complex coordinates 



t = a + iv'v 



(2^ 
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with 
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(29) 



The metric is still defined in terms of the complex 1-form valued in the Jordan algebra 

Dt = p{s)(v'^^^dtv~''\ (30) 

but instead of being given by the 5*1/(3, C*) invariant trace, as 

dsVTrDt*L>t (31) 

with 

't=a-iv^v (32) 

the definition requires to spht C* into R © R as 

Dti + iDt2 (33) 

and the metric is given by the R!^ x R^ invariant trace 

ds^ = Ti Dtii< D\ + Tr £) t^ ^ D't^ + '^dc.+ + "^d^- (34) 

Therefore the metric is not directly determined in terms of the variables t and is not 
invariant under the holomorphic action of the symmetry group of the determinant of the 
Jordan algebra S'L(3, C*). Nevertheless, we will see that these coordinates are particu- 
larly convenient for describing black hole solutions. 

In the same way, the components of C^^'^Eabc: rewritten as a 6 x 15 matrix 

/ P° P3 p2 Q^ p^ Q _p^ Q _^ y Q z q^ \ 

P' Q2 Qs -Qo y P3 ^ q' -z -q" 

2 K. q^ P° P^ P^ Q2 P2 -Pi -y 

g2 -^ -g3 P2 g^ Q^ _Q^ ^ p^ y 

P3 -p2 -z ^ y q^ P^ P^ pi g^ 

\ t P2 z q^ ^y -q^ P^ Qi Q^ -Qo / 

(35) 
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can be recombined into two scalars go = Qo and p° = P° and two Hermitian matrices Q, 
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These coordinates render the A/" = 2 truncation of the theory rather transparent, 
because the latter is simply defined by the real solutions (no £ components) with <;± = 0. 
The coordinates then correspond to the conventional special coordinates associated to 
the special Kahler space Sp{6,'R)/U{3) [41, 42]. Requiring all the symmetric 3 by 3 
matrices to be diagonal corresponds then to the STU truncation. In these notations the 
two truncations of the solutions we will describe in the following will be transparent. The 
corresponding three-dimensional models are defined over the symmetric spaces 

50(4,4)7(50(2,2) X 50(2,2)) c F4(4)/(5L(2) x^, Sp{6,R)) C E,^,)/SpXS,R) (38) 

2.2 Lift to £^7(7) special coordinates 

We will now describe how these coordinates can be generalised to the symmetric space 
E7(7)/5f/,(8) oi^f = 8 supergravity moduh. Just like 5L(3, R) x 5L(3, R) ^ SL{3, €*), 
it is possible to realise -Eece) as 5L(3, O*). Recall that the split octonions are defined as 
doublets of quaternions 

k + is , e + iy (39) 

which are multiplied according to the Cayley rule [46] 

(e + iy){k + is) =ek + sy* + i{e*s + ky) (40) 

and for which the complex conjugation reads 

(k + is)* = k* -is (41) 

For the case of the -E'6(6) truncation one considered the split complex, and these relations 
reduced to say that i"^ = 1 and i* = —i. As opposed to the case of the associative com- 
position algebra, there is no representation of 5L(3, O*) on 3-vectors of split octonions 
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because of the lack of associativity. However, the associativity problem is resolved in the 
case of the representation on Hermitian 3 by 3 matrices over the split octonions. More 
precisely, the action of traceless 3 by 3 matrices over O* on a Hermitian 3 by 3 ma- 
trix closes modulo 02(2) automorphism transformations. This Lie algebra representation 
exponentiates to a group representation of -^6(6) • 

Let us consider a 3 by 3 matrix over any composition algebra which diagonal compo- 
nents are all real 

/ hi 63 e*2 \ 

X= f* h, ei (42) 

V /2 /r h J 

By definition, the infinitesimal transformation of a Hermitian matrix a 

da = Xa+ aX^ (43) 

preserves the condition that a is Hermitian. In order to check that this transformation 
defines a Lie algebra representation, it will be convenient to define the algebra in the 
associated differential graded algebra /\5ll{0*), which can be understood as the complex 
of Grassmann variables valued in the co-algebra on which acts the differential 6 as 



6X' 



-C'jkX'X' 
2 



(44) 



where C^jk are the structure constants. The differential is also defined on vectors asso- 
ciated to the representations of the Lie algebra as in (43). In this formalism, the Jacobi 
identity is then equivalent to the nilpotency of the differential S, as it is in the BRST 
formalism. For anti-commuting parameters, one will use the convention that 

(efT = -re* (45) 

and therefore the identities on the algebra get minus signs associated to the grading. We 
define 

/e3/3*+e^/2+/3e3+/2*e2 \ 



5X = X^- 



\ 








/Jes+ei/i+eJ/a+Zie* 







/2eJ+/i*ei+e2/2*+ei/i*/ 



(46) 



where the second term is required in order to ensure that the diagonal elements remain 
real. One computes that 5"^ does not vanish, but gives 

t3-i(ei,/i)e3 (t2-i(ei,/i)e2)* \ 

5^X = ^x{eiJi)X = ^\ {x^^i{e,J,)h)* ti„,(ei, /i)ei 

t2-i(ei,/i)/2 (ti_i(ei,/i)/i)* 



(47) 
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where the indices of ti{a, h) are defined modulo 3 and 

ri(a, h)a = —-a*{ha) — -b*{aa) 

1 1 

t2(a, b)a = -{aa)b* + -{ab)a* 

Xsia, b)a = -a{b*a + (-l)"a*6) - b{a*a + (-l)"a*a) (48) 

for anticommuting composition algebra elements a and b and where (— )" is minus if a 
is anticommuting {i.e. (— )°' = —1 in (47), whereas (— )" = 1 in (51) to come). These 
transformations generate the automorphisms of the algebra of cyclic products of three 
copies of the composition algebra. Indeed X3{a,b) is manifestly an infinitesimal rotation 
whereas Xi{a,b) and X2{a,b) correspond to spinor transformations. Their group analogue 
are the triality related fundamental representations of Spin{4, 4) satisfying 

{p.{,)l3*a*y = (p.+i(Oa) (a+2(^)/3) (49) 

for the split octonions. We conclude that the algebra closes modulo so (4, 4) transforma- 
tions in this case. This representation corresponds precisely to the graded decomposition 
of e6{6) 

e6(6) ^ 8^-'^ © (8i © 83)'-'' © (gli © gi^ © so(4, 4)) '°' © (81 © 83)'^^ © 8f (50) 

One checks also that £'6(6) is represented in this way on the Jordan algebra such that 

^ X3-i{ei, fi)a3 {x2-i{ei, fi)a2)* ^ 

S'^a=^^ieiJi)a='^ {x3^^{eiJi)a3y Xi^i{eiJi)ai 

' ' \ ^2^t{ei, fi)a2 (ti_i(ei,/i)ai)* / 

(51) 
The coordinates we have been using for SL(Q) will extend straightforwardly to -£7(7) in 
this way. To take care of the associativity issue, we will define the ordered exponential 

: exp[X]aexp[X^ : = a + Xa+ aX^ + - {X{Xa + aX^) + {Xa + aX^)X^) 

+ -(x{X{Xa+aX^) + {Xa+aX^)X^) + {X{Xa+ aX^) + {Xa+ aX^)X^)X^)+. . . 

(52) 

Special coordinates for Et(j)/ SU^{S) 

The GL(1) x£'e(6) coset element will be represented by 16 scalar fields (written collectively 
<;"^) parametrizing the quaternionic symmetric space S'0(4, 4)/(S'0(4) x 5*0(4)), three 
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scalars 0i and three O* valued scalars ^, r/, ( such that the group is represented in the 27 



as 



a = P(^) 



/e-^i \ 
e-^^ 

\ e-^^/ 



o\ 

exp I ^* 
C V* 0/ 



cexp 



r] 
\0 



/e-^i 





o\l 





g02 





V 





6-^3/ 



(53) 



One can still define the complex coordinates t by acting on the identity of the Jordan 
algebra in the conjugate representation 

0\ /e-^"^! \ /0-^-C*\ 

6-2-^2 

e-^-^y 



t=a + i:exp \ -^* 
-C-r]*Oj 



exp 



-r/ 
\0 / 



/ 



a + i 



-201 



-201 



(c-i«.) 



(54) 



\ 



V 



_g-20if* 

-e-2^^(c-ir,*r 



-202 



-201 If 12 



lep 



-202 



-202, 



V + e-^'^'cic-l^v) 

7]* +ie"2'^i(C-|»7*r)? e-2</.3+e-2<^2|^|2+e-2^i|C*-i€»7p / 



The metric is defined in terms of 

/e'^iQ 0\ /000\ 



Dt = p{^) 



e*2 
\ e*^y 



: exp 



TOO 
VC ^* 0/ 



(it exp 



/ ^ C* \ 

r/ 
\0 / 



/e0i 





y 





g02 





V 





£03 y 



(55) 
(56) 

(57) 

(58) 

^''^ 

The coordinates can be chosen independently of a specific parametrization by requiring 
that this determinant is the product of three strictly positive fields. Under this condition 
one can define the Jordan inverse of ft — t) as 



similarly as in (34) by 

ds^ = Tt Dili. D\ + Tt Dtzi^ DI2 + gij{<;)d<;^d<;-^ 
The determinant is defined in the Jordan algebra by [46] 

det[c] = 010203 — %^aj|ajp + aia2«3 + agOgC^^ 

i 

and we have that 

2det[t-t] = 8e-2'^^e-2'^^e-2'^3 

This determinant defines the symmetric cross product such that 



a X a = a — Ti a ■ a Tr 



-Tr^ 



[t - t) ^ = det[t - t] ^(t - t) X {t - t) 



(60) 
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For a Jordan algebra defined over a field, like in A^ = 2 special geometry, the scalar field 
metric would be simply given by [42] 



ds' 



det[t 



:Tr it 



t] X 



[t - t)dt{t - t) X (t - t)dt 



(61) 



which implies that the metric is invariant with respect to a linear action of the symmetry 
group preserving the determinant. However in our case the metric is defined as in (56) and 
is not preserved by the complete Eq^q) symmetry, but only by its SO {A, 4) subgroup. The 
solutions we will describe in this paper satisfy nevertheless that Dtz is real, such that the 
remaining coordinates admit an enhanced linearly realised symmetry 5't/(2) x^2 5'L(3, H), 
as for the solutions sitting in the maximal M = 2 truncation of the theory. ^ 

The electromagnetic charges in the 56 

The electromagnetic charges will be described in the same way. We consider the decom- 
position of eletromagnetic charges according to the -£^6(6) C £'7(7) subgroup, such that 
they decompose into 1 © 27 electric charges, go and the Jordan algebra element 



CL 



( 



\ 



Qi 



+ 



72!±£ 



2 ' "- 2 
Q2 



i'-y 



+ 



9^+1 \ 
■ 2 



79^+i/ 



I 



q^+\. 



Q3 



(62) 



and 1 © 27 magnetic charges, ]P and the Jordan algebra element 



1 



( 



\ 



Pi 



PS-Z _|_ pP3±£ 
2 "^^2 



P2-y 



+ 



9P2±y \ 
- 2 



P3-Z 
2 

P2-y 



9P3±z 
■ 2 



p. 



+ 



?P2 + i/ 



Pl-K. 
2 



I Pl + K, 
' 2 



(63) 



/ 



where g*, Pi, ;\:, z/, z, ^, y, z are quaternions. 

It will often be useful to decompose a 3 by 3 Hermitian matrix over the split octonions 
into a 3 by 3 Hermitian matrix over the quaternions and a 3 vector over the quaternions. 
For this purpose we define the Hodge like operation -k from a 3 vector over the quaternions 



^Note however that this symmetry is incompatible with 5*0(4,4). Indeed, extracting Im[t2] from the 
special coordinates and combining it with the "j's to parametrize the symmetric space i^4(4)/(S'p(l) x 
5^(3)) permits to compensate for the linear £^6(6) realisation, however, only its SU{2) x^a SL{3,M) 
preserves the condition Im[t2] = 0. 
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-1^3 Fi 
V Y2 -Fi / 



(64) 



(65) 



y to a 3 by 3 antisymmetric matrix over the quaternions 

(Y^\ / Y^ -Y2\ 
* Y2 
\YsJ 

Using this notation, we will write 

Cl= q + i-kq , ^ = p + ii^p 

One has then ^ 

(IXCI= qx q + qq^ -Eicqq (67) 

where qq^ is the Hermitian matrix defined from the tensor product of the vector q with 
its Hermitian conjugate g^, and qq is the ordinary action of a Hermitian matrix on a 
vector. In particular, we will consider pqq as the repeated action of the matrix such that 
the factor pq is the ordinary product of two Hermitian matrices, and not the symmetric 
Jordan product. In this way 

det[(l] = det[q] + {q^ qq) 

and 

Tr [Q,lP] = Tr [^p]-(gtp)-(ptg) 

Using these equations, the quartic invariant [47] 



(68) 
(69) 



h = P° det[Q,] - go det[!P] + Tr [(Q, x Q,)(fP x ^)] - - (goP° 



Tr [dT])' (70) 



decomposes as 



h = P%det[q] + (qUq)) - go(det[p] + (p^pp)) + Tr [{q x q){p x p)] 

+ i<i\p X p)<i) + ij>\i X i)v) + {q^v)ij>U) - ij>^pn) - (qUpp) 

- ^(go/ + Tr [qp] - (q^p) - (p^q))' (71) 

Let us mention that all these algebraic relations can be written in terms of SU(2) x 
SU*{6) tensors {SU*{6) = S'L(3,H)). The charge q can be written in term of an an- 
tisymmetric tensor of SU*{Q) Qab, and q as an SU{2) x SU*{6) tensor Q", while fP is 



^Whereas the Jordan product decomposes as 



Cl^ = q^ + qq^ - (g^g)l + H{Tr [q] - q)q 



(66) 
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defined in terms of tensors transforming in the conjugate representations P""^ and P". In 
this way one has 



Tr [d^] = ^Q,,P^' - ]^QIP: 



:(ix(ir=o^ 



_abcdef 



QcdQ 



ef 



e'-^QaQi 



det[Q,] = ^e'''>'<'''fQabQ,dQef + \e''^QabQlQ 



a /^b 



and similarly for P and P'^ 



(72) 



The central charges 

The 28 complex central charges Zij{q,p) decompose accordingly into two real components 

1 



2o 



^/det\vhJ\ 



(qo + Tr [a Q,] + Tr [a x afP] - det[a]/ 



(73) 



and two Jordan algebra elements 

^ydet\vhJ] 
Z^ = p{q) v^det[i;tt;] 

where v is formally defined as 



v[(l+ 2ax T - ax ap°]^^^ 



V 



V 



t-i 



[T- ap^]'. 



(74) 



V 



\ 



fo-^-C\ 



exp 



-T] 
\0 / 



(75) 



/e-^i 

e^^^ 
\ e"'^'J 

with the product appearing in the exponential understood to be ordered such that each 
element acts on the charges in a first place, according to (52). As for A/" = 2, it will be 
more convenient to define the central charges in a complex basis. These combinations 
will break the symmetry group to the one of the maximal M = 2 truncation, i.e. to 
SU{2) xz2 Spin*{l2) C ^7(7) and S{U{2) x f/(6)) C SU{S). Zo and zP are only invariant 
with respect to the subgroup Spi^A) C SU{S), and the traces of Z27 and Z^f 

1 



Tr Z27 = 
TrZ^: 



^det[t;tt)] 

1 
ydet[i/k^ 



Tr [{v^v){(l+2ax T- ax ap^)] 
Tr [{v^v) X {v^v){'P- ap"^)] 



(76) 



20 



are only invariant with respect to the subgroup Sp{l) x 5*^(3) C S{U{2) x U{6)). But 
the 'central charge' 

Z{p, q) = -{Zo + iTi Z^r-Tr Z^ + izP) 

= , .^\ ^ (go + Tr [tq\ + Tr [t X 1 5'1 - det[t]p°) (77) 

is itself invariant with respect to S{U{2) x f/(6)). Of course, if one considers the maximal 
M = 2 truncation of the theory by fixing all split octonions to be quaternions, then this 
'central charge' defines the M = 2 central charge. The unconventional factor of v2 is due 
to the fact that we will not consider integral electromagnetic charges, but rather charges 
valued in a/2Z, in order to simplify expressions. One can similarly recombine the other 
components of the central charges Zjj into the Jordan algebra element 

•DZip, q) = -- (2^7 + iZ2r) + - (z^ + ^Tr Z27 + Tr Z^y - 2Z°)) 1 



P(^) 



^(qot - 2i ■ vQv^ : +iTr [v^v(Z]t + Tr [a(l]l 



-2dei[v'^v] : v^'^'Pv-^ : +Tr {{v'^v) x {v^v)'P]l + Tr [a x a'P]! 
-Ai : v{a x 'P)v^ : +2iTr [{v^'v) x a'P]t 
- ({i det[vh!] + Tr [{v^v) x {v^v)a] + iTi [v^v{ax fl)] + det[c])l 

-2det[i;"f't;] : v^^Uv^^ : -2i : v{a x a)v^ :)/) (78) 

which reduces to the Kahler derivative of the central charge in tangent frame for quater- 
nions. In general it includes 27 complex components, and decomposes into its quater- 
nionic and its i component which transform in the 15 © 2 (g> 6 of S{U{2) x U{6)). In 
particular, its trace appears to be equal to 



-.d 



Tr^Z{p,q) = -^2tdet[t--t]Tr (t - J) - —==== (79) 

L CtJ ,/2?det t- t] 



Z 



It breaks S{U{2) x f/(6)) to its Sp{l) x 5*^(3) subgroup, and it will appear explicitly in 
the definition of the ADM mass of the solutions we will describe in this paper. 

2.3 Extraction of the solutions 

A solution is generated from an Ansatz in the symmetric gauge. For a given nilpotent 
orbit one choses a specific representative h of the semi-simple orbit which characterises 
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the nilpotent element e through [48] 

[h, e] = 2e (80) 

For even orbits, this element also characterises the corresponding nilpotent subalgebra 
as the subalgebra of strictly positive grade, generated by the elements 

[h,e<'''] = 2j>e'''^ (81) 

where < p < n for a finite n. We then decompose a general vector L in this nilpotent 
subalgebra as 

n 

L = Y, L^"' (82) 

p=i 

such that 

n 

[h,L]=2^pL'-) (83) 

p=i 

A solution is obtained from such vector L defined as a function over the punctured R^, 
which punctures correspond to the black hole horizons. The functions defining L satisfy 
a set of differential equations as explained in details in [34]. 

In order to render this section slightly more explicit, let us discuss the particular case 
of the non-BPS composite system in the -£^6(6) truncation. In that case, the relevant 
generator h (which can be expressed in a Cartan basis as [0-0-2+0] [49]) decomposes the 
sp(8,R) algebra into 

5p(8, R) ^ 6'-'^ © (2 ® 3)^-^^ © (g(i © st © s^a) ^"^ © (2 © 3)'^^ © 6^'^ (84) 

where the nilpotent generators satisfy the algebra 

[E:,Ef]=£"^E,, (85) 

The 42 representation of 5*^(8, R) decomposes accordingly into 

42 = 2<-'^ © 3'"" © (2 © 6)'^'> © 8'"> © (2 © 6)*'' © 3^'^ © 2^^' (86) 

where the nilpotent generators satisfy 



I L,,,^ uu,, _ ^ ace '"''■— ■■— -■"'■' 



ab .cdl _ . ^ace^Mf^^^ ^ [j.,,, e^^ = -{S^S^ + SX)ea (87) 
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and the grade 2 generators e*^ commute with the others, but the generators E°. The 
complete C6(6) algebra decomposes as 

^6(6) — 

2<-''©(3^3)(-''©(2®3®3)f-^'©(0[i©sl2©st3©sl3)™©(2®3®3)W©(3®3)('>©2(^> 

(88) 

and the fundamental representation as 

27 ^ (3 ® 1)'-^' © (2 © 1 ® 3)'-^' © (3 © 3)<°' © (2 © 3 © I)*'' © (1 © 3)^'^ (89) 

It follows from the grading that the generators e'^ are nilpotent of order 5 {i.e. (X"^e^*)^ = 
VX^j), the generators e*^ nilpotent of order 3, and the generators Bq, nilpotent of order 
2. Therefore L^ = 0. 

It will turn out to be useful to break the SL{2) symmetry rather drastically by taking 
an Ansatz of the form (breaking also the apparent 5*17(3) to 5*0(3)) 

L = {6ab + Lab)ef + Kabef + Y^e^ 

+ (m - ^S'^'Kab + ^e^^^'e'^fidabLcdKef - L,bLcdKef))ei 

-lil + V- \e-e^^^{23.bK.,K., - L.bK^K.,)y (90) 

The algebra is such that the system of equations can be written for L as 

2 
di<dl. = — [rfL, [L, .^c/L]] (91) 

o 

and some algebra permits to obtain that Lab, Kab and Ya are all harmonic functions, 
whereas V and M satisfy 

di<dV = -e'"'^e'"^^LabdKcd ^ dK^j , di<dM = e'^'^e'"^^ LabdL^d * dK^f (92) 
The Ansatz for the solution is then defined for the E^^Q^/Sp^S, R) representative as 



VtiV*ti 



X e-2^M^^ X 

V X e-^^M^^CoBC X / 

: exp(-2L) = 1 - 2L + 2L2 - -L^ + -V (93) 
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where M^^ = (v^v)^^ . We first focus on tlie matrix element e "^^ M^^ . Computing tlie 
determinant of this matrix gives the scaling factor 



e-i2f^ = det[e-2^M] (94) 

One computes in a specific frame that it is given by the cube of 

^-Au ^ v^-e^'''e"'f K,L,dLef - ^e'^^'^e^'^^L.^L.^nF; (95) 

for the non-BPS solution. It is then straightforward to write an algorithm which extracts 
from M^^ the fields defined in (17). We obtain in this way their expression in the case 
of the non-BPS solution 

t^^ = K,^ + (-M + ie-^^)K, ^ ^^^^^ j£!!!£!j^!Z:^y^ (96) 

where the notation is justified by the property that the real part of t is a symmetric 
matrix, whereas its imaginary part is antisymmetric. The two scalars q± are simply zero 
everywhere. Note that the specific Ansatz (90) has been chosen in order to get such 
simple expressions. In general it is necessary to implement non-linear reparametrizations 
preserving somehow the grading in order to get a solution which is manifestly symmetric 
with respect to the largest available symmetry group. 

Once the components of f* have been extracted, the matrix M is trivial to invert, 
and one can extract the electromagnetic scalars from e~'^^ M^^ Cj:,bc ■ 

It remains then to extract the three-dimensional vector fields. This can be done 
straightforwardly by computing the Noether current 

diy^^VPV-i = -5^-^-^adL'>rfL , (97) 

fc=o '^ "•" ''■ 

In the case of the composite non-BPS solutions 

i,dW = -dL + [L, rfL] - -[L, [L, t/L]] 

o 

= -rfL,,ef - dKa,{ef + ^"^^^'^e.rf - <5'^''(<5'='^e,, - h,)) - dY^e'^ 

+ l{dV+ le-^e"'^/ [Lakd{K^K,f) - Ka,K,,dL,f) ) e^ 

- (dM - ^£'^'=%^'^^L,bL,,dA'e/) ei (98) 
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One obtains 



duj = Ti EdW = ^e"^^ ( da - ^C^^^dCABC ) , (99) 



which in this case gives 



^dcu = dM- ^e'"''e'''fLabL,ddK,f (100) 



and 



dwABC= -^Tr EABcdW 



: e-2^M^^M^^M^^ ^ rfC''''^ - UBcduj , (101) 



which gives 



i.dw"^ = i^duj 5"^ + e'^^^e^'^f {LabdR^d - Kabdl^d) + le^^^dY^ 

-kdVab = dLab 

i,dv^ = -i,duj (102) 

In order to restore the manifest covariance, we moreover carry out gauge transformations 
hnear in time in order to remove the -kdu components as well as the corresponding non- 
covariant constant terms in 

C° = 1 - e^^^e'^^^e^'^^L.feL.rfLe/ (103) 

such that 

C^dt + u) + v^-dt=(l- e^^-e'^^^e^'^fLabL^dLef) (dt + cu) - (dt + cu) 

6 

= _e4^l£»-£^'^/L,,L,,Le/(rft + u) (104) 

and similarly for C"^(dt + uj) + w"^ + did^H). 

Once the fields of the truncated theory have all been extracted, one can simply re- 
organise them in the Jordan algebra scalar field t and equivalently for the vector fields. 
It appears that all the nilpotent orbits of e6(6) admit at least two grade zero components 
associated to the two extra generators of the Cartan subalgebra of e6(6) with respect to 
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the Cartan subalgebra of 5p(8,R). This suggests that one could chose the coordinates 
such that the two scalars <;"-!- = in the Ansatz. Equivalently for A/" = 8 supergrav- 
ity, the maximal nilpotent orbit of -^8(8) that admits a non-trivial intersection with the 
coset component has a grade zero component which is defined by the non-compact gen- 
erators of the algebra so (4, 4) C e8(8). If the two scalars ^-t = in the truncation, 
it seems reasonable to believe that the 16 scalars ^^ that parametrize the quaternionic 
space S'0(4,4)/(S'0(4) x S'0(4)) can also all be set to zero. However the situation is 
not always so simple, and preliminary computations show that the scalar fields <;± get 
non-trivial values in the maximal solvable systems. It seems nevertheless that the scalar 
fields parametrizing the submanifold G2{2)/SO{4) are still constant, such that one could 
obtain the N' = 8 scalar fields in S'0(4, 4)/(S'0(4) x S'0(4)) by promoting the two non- 
trivial scalars e^'^='= to uni-modular split octonions. Anyway, the systems we will study 
in this paper have more symmetries, and in these cases the two scalars i;± are indeed 
null. Moreover, one can then check that the scalar fields q^ are not sourced by the other 
fields, through a cancellation of the contributions of the other scalar fields and the elec- 
tromagnetic fields. We will call the scalar fields that are constant in the solutions the 
'flat directions', referring to the property that they decouple completely from the equa- 
tions and are simply set to constant. Note however that they do not define the common 
stabiliser subgroup of all the electromagnetic charges of the black hole constituents. 

In all the systems there is a linearly realised [5^/(2)]*^ symmetry. This symmetry per- 
mits to determine the ordering ambiguities appearing in the generalisation of an expres- 
sion written in terms of real functions to the expression written in terms of quaternions. 
Each quaternion transforms with respect to one specific SU{2) on its left, and one specific 
SU{2) on its right. The only problem that can appear, is when quantities like 

xy* — yx* (105) 

appear explicitly in the solution (where x and y transforms identically with respect to 
[SU{2)]'^), since such expression vanishes identically for reals. In fact one might have 
to deal with such ambiguities for the most general systems, and it may then become 
necessary to compute the solutions within the Ej(^t^ / SU^{4, 4) non-linear sigma model. 
However, we will mainly focus in this paper on solutions admitting a much larger linearly 
realised symmetry SU{2) x SL{3,M), such that quaternions all combine into 3 by 3 
Hermitian matrices, or 3-vectors over the quaternions. In this case we will see that there 
is no such ambiguity problem. 
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3 Composite non-BPS system 

The nilpotent orbit in which is parametrized the composite non-BPS solution of the STU 
model [34] can be embedded in a larger c^is) nilpotent orbit. The nilpotent algebra admits 
as an automorphism SL{2) x SU(2) x SL{3,Ii). The parametrization of the solution 
breaks SL(2), which is related to Ehlers symmetry. But the automorphism symmetry of 
the Jordan algebra is preserved because the solution can be explicitly written in terms 
of scalars, Jordan algebra elements, and 3 vectors over the quaternions. The nilpotent 
algebra itself is entirely determined by the gii weights and the Jordan algebra cross 
product, such that the satisfaction of the equations of motion in the -Eece) truncated 
theory ensures that they are satisfied in A/" = 8 supergravity. The structure is as such 
extremely constrained, and one can check at each step that the SL{3, R) invariant terms 
generalise in a unique way to the SL{3, H) invariant expressions. 

3.1 The solution Ansatz 

The solution is defined in terms of harmonic functions organised in two 3 by 3 Hermitian 
matrices over the quaternions, /C and C, and a 3 vector over the quaternions Y. There 
are moreover two sourced scalar functions V and M which solve the equations 

di. {dV + Tt Cd{IC X IC) - Tt IC X ICdC) =0 

di.{dM-TT Cx CdlC)=0 (106) 

The system is invariant with respect to the SL{3, H) symmetry group of the cross product. 
The four-dimensional metric is determined by its scaling factor 



e 



-4U 



Vdet[C] - {Y^C xCY)-M^ (107) 



and the Kaluza-Klein vector 

i.duj = dM -Ti Cx C die (108) 

The scaling factor can be seen to be invariant with respect to 5*17(4, H), by rewriting its 
two first terms as the determinant of a 4 by 4 Hermitian matrix over the quaternions. 
This is a consequence of the property that the system reduces for /C = to the sys- 
tem describing single centre non-BPS black holes which is invariant with respect to this 
symmetry. 
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For the sake of clarity, let us recall the form of these equations in terms of SU(2) x 
SU*{6) tensors, as discussed in equations (72), 

d^(dM- ^s'^'^^'^'fLabLcddKef) = (109) 



and 



abode f T T T 

— fc IjabJ^cdJ^ef ~^ 'Ta^'^l^^ ^ a. ^bc^de-^ f 



e'*" = VT^c""'"L,iL,jL., - —Ccfe""^'%'U,LjeYf - hf 



,,, = ,M-ie--/L..L..dAV (110) 

Nevertheless, we shall only use the Jordan algebra notations in the following. 
The expression of the scalar t is 

and the 16 scalars q^ are constant. Here one sees that setting Y to zero, the scalar fields 
are Hermitian matrices over the quaternions and the solution sits in the M = 2 truncation 
of A/" = 8 supergravity. The product with Y is defined as for an ordinary 3 by 3 matrix 

det[£] ^ ' 

The three dimensional vector components of the vector fields are 

* dwo = dV -Ti Cd{ICxlC) + Ti ICxICdC 
*duj = 2CxdlC-2K,xdC + HdY 
■kdv = dC 
i.dv° = (113) 

and they combine with the scalar components to give the 28 vectors 

A"^ = -e^^ det[C]{dt + oo) + v'' 
A = -e*^ (det [C]JC-MC + i^{Cx C)Y) {dt + u) + v (114) 

and their 28 duals 

Ao = e^^{MV + MTr [(/C x /C)£] - l^Tr [(£ x £)/C] 

- det[£] det[/C] + 2{Y^{IC x C)Y)) {dt + io)+wo (115) 
A = e^^{VC xC- 2M/C xC + det[C]IC x /C 

-2C X YY^ + £*(-/C(£ X £) + M)Y) {dt + co) + ^u 
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They can be rewritten conveniently in terms of the scalar fields and M as 



A" = -4 " ^ {^f. + uj) 

^2zdet[t - t] 

A/2zdet[t - tj 

^ = 4^^ _ (t X t)((it + cj) + zv 

^J2idet[t - t] 

/ e^ 

Ao= ^ (Hpt[f] + det[t] - Tr [t X tt] - Tr [t x tt 



>/2idet[t - t] 

v/2i det[t - t] m\ {dt + u) + wq (116) 



g3^ 



3.2 Solving the differential equations 

Let us now solve this system. We define the harmonic functions 

^^^^^ ry rf ^^^^^ ry ry ^^^^^ ly ry 

Note that we already fixed each pole of /C to be proportional to the corresponding pole 
of £, because this is required for the Noether charge associated to each centre to be 
nilpotent of order 3, such that it can correspond to a regular black hole solution [34]. 
Although, the solutions for V and M can be found in principle, they do not admit closed 
form expressions in general. Therefore we will first compute these functions locally near 
each pole Xa, and in the asymptotic region. We will discuss the global solution when we 
will compute the angular momentum. 

Near the pole x^, C admits the expansion 



rr* \ ^^^^^ nr* ry I \ ^^^^^ ry ry O / 

I ■^^l njzA I -^ Xa\' B^i/l I -B -^1 

+ — I > - -=■ \6\Xb X a)%\^B '^Aji Oj) ps ^ a\ ) \\^ '^a) \^ ^ a) 

2V^^^ \Xj.—Xa\ 'J 

+ 0{{x-XAf) (118) 
where all higher order corrections are symmetric traceless polynomials in x — x^. /C 
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admits the similar expansion 

/C = 7a 



PA ,^^+^- ?^ 






7bi 



+ ^(E^^i 



Bt^A 



•Xj ri tXj I 



+(5: 



7b 1 



PB 



Bt^a 



X R X I 



-\Xq Xj^ji 1 (^X Xj^j 



PB 



2\ ^^^^ nrt nrt O 

N \JU'R ^ A \ 



(^d(Xs XAJiiXs ^AJj Oij\XB Xa\ ) j [^ ^AJ [X ^a) 



+ 0{{x-x,^f) (119) 
Using this expression, we find the particular solution for M near x^ 



M^^' = 7a ■ 



det[p/ 



\x — X, 



+ 7ATr \pA ^PaU+Y^ 



PB 



Bf^A 



\X fj X A 



\X X I 



+ ( E ^-Tr [p. X p, p,] ^^^3^4^) Ir^ + 0(1) (120) 



Bt^A 



<T^ O / /y /y 



to which we can add a homogenous solution 



M(°' 



m 



P\ 



E PA _ sr^ 
'f -y ^ -^ 



J j\\X Xj^Ji 



'1211 



rtr> ^^^^^ rf> /y 

A ' -^l A ' -*' 



No higher order poles are allowed because they would necessarily produce singularities 
in the metric. The absence of Dirac-Misner string singularities requires to constraint the 
simple poles of M'°* such that 



Pa = 7ATr 



PA 



fH-E 



PB 



Bf^A 



\X u X / 



'^Y. 



PC 



C^A 



\Xc Xa\ 



and 



^ — £ijk 



J\{x — XA^dx^ 



lAj lAj A 

One finds therefore the expansion of M near Xa 



+ 0{x - Xj 



(122) 



(123) 



M = 7A^-^^ + 7ATr 



\X "^ a\ 

Jj^[X — Xj^ji 
I I '\ 

\Jy ^ A \ 



PA X paU + 22 



PB 



+ 7ATr 



PA 



'^E 



B=/=A 
PB 






\Xp! X I 



X 



' + E 



PC 



B^A 

+ ( ^ 7i?Tr [pa X PA pb] 



IXn X / 



C^A 

,Xg ^AJi \ [•^ ■^a) 



\x — X, 



Bf^A 



I r^ / I 19 



C(l) (124) 
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We will solve similarly V by rewriting its equation 

di<d{V -Ti CIC X /C) = -2rfTr /C x /C i< dC 
We find in the neighbourhood of x^ the particular solution 



(125) 



(y-Tr /:/Cx/C)'" = -2 7^ 



det[pA] 

I 1 3 

Jj ^ A 



+7ATr 



Pb 



PA^PA[k.+ ^IBT—^ 



B=/=A 



B j4 



+ 7i ( E Tr [pax PA p^] 1^^-^)1^-^ + 0(1)1 (126) 



B^A 



to which we can add the homogenous solution 

(V-TvCKxICf'^h + Y, *" 



E 



/J^(X X^jj 






(127) 



The definition of the associated charge Qqa through the explicit form of the vector wq 
requires that 



Qqa = qoA - 2Tr 



PAyl<i+ ^IB- 



PB 



Bf^A 



ih f> Jb I 



X 



[k,+ ^lc- 



PC 



C^A 



Xc Xa\ 



such that 



, , goA , ,/ P\{x-XAydx^ 
dwQ = -kd- + d I Eijk j ^ h 0[x — Xa, 



■X> iX-' A 

The expansion of V near Xa gives finally 



V = ll 



det [p/ 



\x — X, 



Tr 



7JTr pa^PaU +^ 



2^apa[[+J2 



X — X, 



PB 



:i28) 



(129) 



Bi^A 

l3\{x-XA)i , goA 



rtf> I rf /y» 



Bi^A 



\% '^ A 

PB 

Jb R •Xj I 



\X 'C A 



X 



V r-< \Xr — Xa\^ 



-pA 



t+E 



7b 1 



PB 



Bi^A 



\Ju^ 'J^ a\ 



Cj^A 



^+E 

Ct^A 



Ici 



PC 



M-'C "^ A 



\x — X, 






( E 7a(27b - 7A)Tr [p^ x ;.aPb] j^^^l j^^ + 0(1) (130) 



^B^A 



I I "? I I 19 

S A / A 
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3.3 Near horizon geometry 

Let us now study the form of the scahng factor e~^^ near Xa- We get 

V^det[£] -M^ = ^^M^(/3^ + 2^aJa)\x - x,), + 0{\x - x,|-^) (131) 

X ^ A 

Therefore regularity requires 

f5\ = -2^aJ\ (132) 

Using this equation one then gets 

\X — Xa\ \ 



Tr 



\J\{x-XA)i) ^ ^^^1^ ^ |_3 



•^ A 



0{\x-Xa\-^) (133) 



To interpret this hmit we must compute the charges ({a associated to the vector field u). 
Near x^ one computes that 

CxdlC-lCxdC= {-lAt - t+ y (7A - 1b)-, — — — r) x rf , ^^ , + Oil) (134) 

and therefore 

qA = 2pA X (iaI: -H+y^ilA- 7b) I ^"^ 1 ) (135) 

We then use the Jordan algebra identity [46] 

det[pA]?A = 4pA X {qA x {pA x pa)) - pa x P^Tr [;7a?a] (136) 

to show that 

-qoAdet{pA] + Tr [{qA x ^a)(pa x pa)] - ^(Tr [pa?a])^ 



Fdet[£]-M' 



\2C '^ A 



^^^^r^^-^^ + Oilx-x^n (137) 
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We recognise the first term as a component of tlie quartic invariant. Using coordinates 
such that J\di = Ja^z "we can rewrite the scaling factor as 



g-4C/^ 



-qoA det[p^] + Tr [{qA x qA){pA x pa)] + {q\ (pA x pa) Qa) - j{^^ [pa^a])' 



X '^ A 



,2 



^■^^^"^ '''^'^- + 0{\x-xA~') 



CC ^ A 



. '^"y^^^P^^ + JA^o^y^.^ ^ ^(1^ _ ^^|-3) (138) 



\x — X 



A\ 



as expected for a non-BPS extremal black hole [40, 50, 51]. The horizon is therefore 
squashed, and the scalar fields do not take constant values, but rather 

*^ = ^^ ^^^ ^ ^W. (139) 



dqAO 

This corroborates the analysis in [32]. We will see in the next section that the scalar fields 
satisfy a generalised attractor equation which is the U{1) Ehlers transform of the standard 
one, for which the angle depends on the angular momentum J a and the coordinate 9^- 

3.4 The asymptotic region 

Let us now consider the solution in the asymptotic region. In that case one has 



|a;| yxY \x\ \xy 

(140) 

The boundary conditions for the functions V and M were fixed at the poles, and without 

solving the equation completely one only fixes these functions up to arbitrary harmonic 

functions. Nevertheless, the conservation of charges determines the first pole to be the 

sum of the poles associated to each black holes and we can determine that 

A# , Tr [r X lY^AlAPJ^ , ^^,1 2n 

M=mA j— -^(D{\x\ ) 

\x\ 

y = /z+^-^^°^ + ^^[''"^^^|^^^^-^"^^^^-^]+0(|a:n (141) 

One can check that this is indeed the case for the global solution we will discuss latter. 
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The ADM mass 

Using these formulas one obtains that 



-4U 



hdet[[] - {y\C x [)y) - m^ + (hTr \[ x [^pA 

^ A 

+ det[/:] (^J2 ^OA + Tr [2/: X liJ^^^P^ - k.x C^Va 

A A A 

- {y\C X C)Y,u) - ((E^-4) V X t)y)-2{y^[[ x ^p^)i/ 



— 2mTr 



[xCY,1aP^\t^,+0{\x\-^) (142) 



The conventional boundary condition on the scaling factor of the metric determines 

1+ [y\[ X [)y) +m^ 



h 



det[r] 



To interpret this formula more precisely we will also need the asymptotic moduli 
and the total electromagnetic charges 

% = ^%A p = ^PA q = ^qA 

AAA 
A A 

Let us first combine the terms 



(143) 



(144) 



(145) 



ra 



det[r] 



Tr 



Cx[^pA+ det[r] {V qoA + Tr \2[ x ^ J^TaFa - ^ x I^Y, 



.PA 



-2?TiTr 



i'-'E 



lAPA 



:det[r]Tr 



(^-d4l"^)"(^ 



det[r] 



ra 



E 



PA 



Aei[C]Ti 



L J ^ A 

\ det[r] / V^ det[r] r V^ dei[[] /\ 



(146) 
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We can now have a look at the terms in y 

^^^^^^j^Tr [CxCp]- {y\[ x [)q) - {q\[ x [)y) - 2{y\[ x j,)y) (147) 

For the Jordan algebra associated to the associative composition algebras, we have the 
matrix product identity 



p{C X [) + 2C{[ X p) = Tr [r X r p] 1 



(148) 



which is generally valid for any exceptional Jordan algebra for the commutative Jordan 
product. One can then simplify (147) as 



det[C] 



[ X [ \ t [ X [ 



y P: 



-\y 



[ X f \^ 

■y] Q 



[ X [ 



q' 



-\y 



Using the octonionic expression for the asymptotic moduli and the electric charges 



(149) 






detfr] 
we can combine these expressions in 



det[C]' 



0,= i + i-kq 



(150) 



Tr 

and 



[(^- 


[ 

det 


r 




Tr 


[( 



m X U, 



[ 



detf/"] 



m\]} 



k.- 



l 



-m) q 



[ X [ \\ 

y] q 



[ X [ \t /■ X /■ 
[x [ 



det[[] )^ VVdet[/'] 
This way we can rewrite the scaling factor as 



-w 



l + detlrlUo + Tr [floQ,]+Tr 



I 



q' 



I 



det[r] 



y 



:Tr [floXflop] (151) 
Tr [floQ.] (152) 



det[/'] det[r] 



+ flo X flo P 



— \0{\x\ 
\x\ 



(153) 

We conclude that the ADM mass is determined in terms of the asymptotic moduli and 
the total electromagnetic charges as 



M. 



det[/'] 



go + Tr [a„Q,]+Tr 



I 



X 



I 



LVdet[/'] det[r] 



+ flo X flj, p 



This can be written in terms of the fo themselves as 



M. 



1 



go + -Tr [(io + f(,)Q,] + Tr [i,, ^ t^V 



(154) 



(155) 



v/2zdet[/o- Y\ ^"" ' 2 
Although this solution does not sit within an A/" = 2 truncation, it is characterised in 
this duality frame by the 'central charge' (77) such that 



M^oM = ^(-Z(g,p) + Tr [®Z(p,g)] 



(156) 
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The angular momentum 

This perturbative expansion does not permit to compute the angular momentum. In 
order to obtain the global solutions, one needs to solve several Laplace equations with 
sources. The ones involving only two centres were already solved in [14], and we will 
refer to this paper for further details. However, the equation involving three centres 

AFa,bc = I ^ | A (- 1 -) (157) 

\Uj "-" A.\ \ "^B"^ "-" C \ / 

was only solved in the axisymmetric case for which Xa,Xb, Xc were assumed to be aligned. 
In this case 



F 






A fif / \ / Nil II II II II II 

■^1 ^^ ( ryt ryt i , i /y 'IT' \ \ /y /y /y /y /y /y /y /y /y /y /y /y 

KJUq "^ a. } \'^ C ^ a. ) \ \^ ^ B \\"^ "-" C \ l"-" B ^ C \ \ "-" B \ \^ B "^C" '^O 

(158) 
We will not be able to solve this Laplace equation in general, but one can define 

Fa,bc = - / 1^ I n n m \3 ^^^ 

J 2n \y- x\\y- XaWv- XbI-^IV - Xcl'^ 

which is the unique solution to the Laplace equation that admits no poles and is regular 
in the asymptotic region. We study this function in the Appendix A, where we prove that 
the integral converges, and we calculate its Taylor expansion in the asymptotic region. 
The exact solution for the function M is then 

, V^ / det[pA] , Tr [ PA X PA , Tr pA^ x [\ ^^^ JAi{x - x^Y 
M = m + > 7a I ^ + —. ^— + -. r- - > 



/y I O I ly /y I ^ I /y ly I / ^^^^^ I /y» /y» O 

^ X, A| \ A\ I A| / ^ I A| 

+ y"7ATr [pAXpni-. n T + i m T " i m 1) 

^^^^^ \ ly ly ly ly rf ly ly ly rf rf rf rf I 

ViX/ '-*^A'-*^ ^ B\ \"-' a ^ B\\^ ^ A\ \^ a "^ B \\^ "^ B \ / 

Iv^rp ^ f lA + lB , lA-lB (\x-Xb\ 1 



2/ J L ^ £ J^ L ^~^ \ \ |'>| II I9\l 19 I 
^^^^^ \ rf rf i^ rf rf rf rf \ i^ \ \ rf rf i^ rf rf 

\J/ "^A"^ "^ B\ \"^ a "^ B\ ^i-tj "-'A \^ Jo j^ 

+ Vl 7cTr pAPB X pc'(-^A,BC + 1 n n r) (160) 

^^^^^ \ rf rf rf rf rf rf I 

N \Jij A Ju r^WJu -D Ju c WJ-' J-' C \ ' 

A^B^C ^ ^ cii B CI I CI 

We will not compute the Kaluza-Klein vector u globally, but the asymptotic behaviour 
of Fa,bc at X — !■ oo is enough to compute the asymptotic behaviour of the vector 

-kduABfi = d[Fi^A,B)c + 1 n n f ) ~ 1 n r^i r (^^i) 

y ^ ' rf rf rf rf rf rf / rf rf rf rf rf rf 

^ \Jj A Jj f'' \ \Jj fi Jj (~- \ \Jj Jj (-• \ ^ \Jj Jj A \ \Jj J/O \Jj Jb (-• \ 
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as 

1/ /y>^ rp^ rp'' /y>* /y> ^ I rp'' / T* \ tT . . , rpj/'frp'^ 

^AB,c = o 1 n 1 + I n r + 1 n r — m — 



12 / i 

rp rp I ■" I /T'l'' 



(/Tri^ rp '' \ ( rp rp | . ( 'Y' 'T* I I 'T' 'T' 11^ i 'T't' /i '~f 

B "^cJ y-^A -^Bj y-^ B ■'^Cjy.'^A "^bJ I'^tlk-'^ ""^ 



<T^ O /T^ <T^ /T^ ri/-> ry ry I ry ry I ry ry I /t^ ry \ 



(162) 



The solution is 

J\{x — x^ydx'' 



UJ 



^^^'^ 1^_^ 13 ^ Z^ ^C-T^ /'^PS ^ PC ^AB,C 

ilA - 7B)Tr /■ PA X pB eijk{x\ - x%){x^ - xi)dx 



A ' ^' A^B^C 



^— ^^ /Y« fY» (^ fY» ry ry I /y* ry I ry ry I /y* /y* 1 

|X^ JyB||X X_4||X Xa|^^|X X^l -r |X Xb| -r \JyA -^B]) 

sr^, _ , Tr pAVA x pBSjjkjXg - x^^)(x-^ - xi)(jx^ 

/ V V '^ /-Bj i TTi Tol i 

^^^^^ rp rp i( fv* rp i^ rp rp 

Using the asymptotic expression of lOABfi one obtains the angular momentum 

A zl\ R \ C^A 

-(7b^+ ^(70-75)1 ^^ — r) xpBPA I 



(163) 



(164) 



^>B ^ Ct^A 



C^B .-. ^^.,. , :Xa Xf 



ry " I ry'' ry '' \ ( ry ry \ , { 1p 1p \ ( ry'' ry <' 1 

+ 2^ 7c -Ir PAPB X PC 



rift rfi ry ry ry rp I rp rp I rp rp I rp rp 

A-J-D-J-C I ^ -^bW-^a -^cW-^b •^C\\\-^A -^B} t \J^A -^Cl ' \-^b -^c 



^^ + 2 5^ ('^°^^^ ^ '^'^ f^^^""^ " '^'^ f^^^""^ " ^^^°'') ■ 



A ^ A>B^ ^1^^ ^^1 

ry " I ry'' ry <' 1 ( ry ry \ , ( OT' Hf | I O^ Ip \ 

+ 2^ -fc i-^ PAPB y< PC-. r. r. rr. r—. r—. ^ 

^^•^^m ry ry ry ry ry ry I ry ry _\_ ry ry _l_ ry ry I 

AfriiQ \-'-'A Xs||X^ J-,Q\\djg •>^c\\\-'^A -^bI ~r l-^A -^d ~r |-^B •'-'C\j 

The first term is simply the sum of the intrinsic angular momenta of the black holes, 
and the second is the standard angular momentum resulting from the non-commuting 
dyon charges. The third one is rather non-standard, and is generated by interactions 
between three non-aligned black holes. One can easily compute that the sum over the 
permutations of ABC of the function 

\X A Xb\ K^C ^ b) \X a ^Bj'y^c Xb)\Xa -^ b) (165) 
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vanishes. Therefore the associated contribution to the angular momentum clearly van- 
ishes if 7a = 7_B = 7c- The function of the charges {jc + c)Tr paPb x PC: does not 
obviously look like it is -^7(7) invariant for any given constant c. Nonetheless we will show 
in the next section that it is. 

4 Duality invar iance 

The composite non-BPS nilpotent orbit is characterised by a semi-simple element h 
of so* (16) (which can be expressed in a Cartan basis as [.oSoooooj) defining the graded 
decomposition 

so*(16) ^ 15'-^' © (2 ® 2 ® 6)<-'' © (g^i © sh © su(2) © su*(6)) "" © (2 © 2 ® 6)'^^ © T5*" 

(166) 
such that 

128 ^ 2^-^' © (2 © 6)'-'^ © (2 ©T5)'-'' © (2 © 20)'"' © (2 © 15)'^' © (2 © 6)'" © 2'^' (167) 

A generic element of the grade one component is not left invariant by any generator 
of so* (16) of negative grade. To compute the isotropy subalgebra in the grade zero 
component, it is convenient to decompose 

su*(6) ^ (2 © 4) <-=''© (0 ti © su(2) © su* (4) )<"'© (2 ©4)'^' (168) 

such that 

15^ 1'-"'© (2 ©4)'-^' ©6'" (169) 

A generic element of the (2 © 15)*^' is then chosen to be a doublet of non-null orthogonal 

vectors of 50(1,5) = S'f/*(4)/Z2 with a non-trivial component in the grade —4 singlet. 

If one vector is time-like, the isotropy subgroup is 

3 
(J[ SU{2)a^ k (0 2= © 2a © r) C SpmlilQ) (170) 

A=0 S>A 

which defines a Lagrangian submanifold of dimension 83 inside the £"8(8) nilpotent orbit 
(associated to the weighted Dynkin diagram [.oiooooo 1) of stabiliser 

Sptn{4, 4) K E6X8+6 c ^8(8) (171) 

If the two vectors are space-like, the isotropy subgroup is of the form 

(5f/(2)i X SU{2)2 X SL{2, C)) x (2^ © (2i © 22) © 4 © r) C SpinKlG) (172) 
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which defines a Lagrangian subinanifold of dimension 83 inside the -^8(8) nilpotent orbit 
(associated to the weighted Dynkin diagram \.^°^^^^^ 1) of stabihser 



Spin{5, 3) X R^^«+^ C ^§(8) (173) 

For single-centre extremal black holes, the regularity of the solutions requires the two 
so(16, C) weighted Dynkin diagrams, determining respectively in which Spin^{16, C) orbit 
and in which -^8(8) orbit the Noether charge sits in, to be identical [26]. Only the first 
representative (with one time-like vector) satisfies to this criterion. In the neighbourhood 
of a given centre, the semi-simple component of the stabiliser subgroup of the momentum 
P G 128 must necessarily be compact in order to interpolate smoothly the stabiliser 
subgroup 5*^(4) k R^^ associated to a regular solution. Therefore this criterion must still 
be satisfied in the near horizon regions of the solution. But in principle the signature 
of the isotropy subgroup of the momentum P G 128 is not necessarily constant allover 
space, '^ and we could have regular solutions for which the total Noether charge defines a 
representative of the other nilpotent orbit. 

A configuration of charges which gives rise to a regular black hole must lie in a smaller 
orbit. After some computation, one finds that the only configurations corresponding to 
regular extremal black holes in this system are necessarily corresponding to non-BPS 
black holes. To illustrate this, let us consider a particularly simple configuration, with 
one antisymmetric tensor in the grade 1 component P^ of weight 1 with respect to a given 
Cartan generator of £[2, and a singlet Qq in the grade 3 component and of weight —1 
with respect to 5(2- Then one computes that the stabiliser of such charge configuration 
only allows for two components of grade 1 and —1 satisfying 

^^"'^aU ^a+"Qo = (174) 

the generators of su(2), and the generators of su*(6) that leave P^ invariant, as well 
as the combination of the weight 2 3(2 generator E'^'^ and the grade 2 component -E"'' 
satisfying 

E++Q0 + P"^P+ = (175) 

and the grade 1 component of SI2 weight 1 E'^""'. In order for the grade component of 
the stabiliser to be compact, we see that P^ must have eigen values of the same sign, such 



*The composite BPS solutions which total electromagnetic charges have a negative quartic invariant, 
admit for instance a mometum P of stabiliser {SU{2) x SU{6)) k C^^^ ® R in the near horizon regions, 
and of stabiliser {SU{2) x SU*{6)) k C^''^ e K, in the asymptotic region [34]. 
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that it admits 5*^(3) as an isotropy subgroup of SU*{6). In order for the combination of 
E^"'^ and -F+" to define a compact generator, we see that we must have 

which is equivalent as to have a negative quartic invariant. The components E'^'^, E""^ 
and E~^°'°-, then combine into nilpotent generators transforming in the 27 of 5^(4), and 
this reproduces the correct isotropy subgroup of a regular non-BPS black hole [37] 

Sp^iA) X R2^ C SpinlilQ) (177) 

The elements of the 2 ® 20 grade component in (167) correspond by definition to 
fields which are not sourced by the system of equations. With respect to 5^(1) x Sp{3) 
they decompose into 2 (g) 14 non-compact generators and 2 (g) 6 compact ones. The non- 
compact generators are associated to the scalar fields, and the compact generators to the 
electromagnetic fields. The 'fiat directions' therefore define the quaternionic symmetric 
space 

F4(4)/(^p(l) x^, Sp{3)) C i?7(7)/^f/c(8) (178) 

Indeed, the imaginary part of the scalar field t is quaternionic in the Ansatz (does not 
include components in i). The corresponding non-compact generators in e6(6) are realised 
as off-diagonal anti-Hermtian 3 by 3 matrices over the octonions which are linear in i. 
The group of unitary matrices over the split octonions acting on traceless Hermitian 
3 by 3 matrices extended by the action of the authomorphism group 6*2(2) defines the 
fundamental representation of -^4(4). Therefore we see that the scalar fields which are 
trivial in the solution indeed parametrize the quaternionic manifold (178). 

4.1 Generalised fake superpotential 

The isotropy subgroup of SU{8) of the element h defining the nilpotent subalgebra is 
^^(l) X 5^(3) (because the U{1) factor lies in the Ehlers group). As it was discussed in 
[34], this implies that the set of inequivalent embeddings for the solvable system described 
by the 44 functions C,IC,Y,V,M is parametrized by the 39 angles parametrizing 

SU{8)/{Sp{l) X Sp{3)) (179) 

which indeed add up to the 83 dimensions of the orbit. These angles parametrize two 
orthogonal antisymmetric tensors of SU{8) satisfying 

Uiju'^ = 2 coiijUki] = UikU^'' + Qikn^'' = S^ oo'^^jk = 

l_e^jklmnp,^^^^^^^^^^^^^ = 1 (180) 
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where u^^ and Q^^ are the complex conjugates of Uij and Qij. In the solutions described 
in the last section, Uij defines the 'central charge' 

u'^Zij = Z (181) 

and flij the Jordan algbera identity, such that 

nijZ'^ = Tr [(DZ] (182) 

The mass formula can be rewritten as 

M^^^ = ^{n,jZ'^-u'^Zi,) (183) 

Note that within our Ansatz the potential QijZ^^ — u^^Zij is automatically real, whereas 
its imaginary part was interpreted as the NUT charge in [34]. This discrepancy is due 
to the fact that we already assumed that the NUT charge vanished when we computed 
the electromagnetic charges, whereas one would need in general to consider the complete 
spatial components of the field strengths. 

Let us define the generalised fake superpotential 

W^^{n'^Zij-uj,jZ'^) (184) 

which within our specific Ansatz is defined as 

W = -{Tt['M]-Z) (185) 

The real part of CDZ is a Hermitian matrix over the quaternions {i.e. the components 
in i all vanish), and rewriting this property in terms of the tensors cuij, Qij we conclude 
that in general 

Uik^ji {Z''' + u'^^n'^Zp,) = (186) 

These 12 constraints (in the 2(8)6 of Sp{l) x 5*^(3)) reduce the number of angles associated 
to a duality frame for given asymptotic central charges to 27. One of these angles is 
moreover determined such that 

^{n,jZ'^-uj'^Zi,)eRl (187) 

In this way the 56 remaining parameters defining the specific nilpotent element associated 
to the asymptotic momenta V~^(iV|i28 are the 28 complex central charges. 
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Relation to the standard fake superpotential 

This potential W somehow generalises the 'flat direction dependent' fake superpotential 
described in [27]. To see this, let us consider for a moment the exceptional M = 2 super- 
gravity with moduli parametrizing the exceptional special Kahler space -E'7(_25)/(f/(l) x 
-^6(-78)) [41, 42]. In this case our Ansatz applies as well, with all Y set to zero and £ 
and K, understood to be three by three Hermitian matrices over the octonions. In this 
case the constraint (186) disappears and 

W = ]-{Ti [fi®Z] - det[(]]Z) (188) 

where 1] is a complex Hermitian matrix over the octonions ^ satisfying 

1] X 1] = det[fi] n , Tr [fiH] = 3 (189) 

VL is parametrized by U{1) x _E'6{-78)/-^4{-52)- And so with the additional positivity con- 
straint VT > 0, the remaining free parameters lie in -E6(-78)/-^4(-52)5 which deflnes a 
compact version of the space of flat directions -E'6(-26)/-^4(-52) of non-BPS charges in the 
theory [52]. Therefore ^^(fi) can be interpreted as a 'flat direction dependent' fake super- 
potential as in [27]. In order to be the case, one should recover the single-centre non-BPS 
fake superpotential when the auxiliary parameters extremize the potential. However the 
constraint that VT > is relatively complicated to enforce, and in order to simplify the 
problem we will consider instead the potential 

w. . "^^g^ (190) 

cos(2aj 

as a function of fi G f/(l) x -E6(-78)/-^4(-52)5 which is equal to W when W is real. It turns 
out that this potential coincides with W at its extremum, if and only if a is the phase 
appearing in the deflnition of the non-standard diagonal form of the central charge and 
its Kahler derivative [26], and in this way, there is no ambiguity in the deflnition of Wa- 
Using the property that the constraint (189) is invariant with respect to the variation 
of VL in f/(l) X £'6(-78)/-^4(-52), one relates the variation of Vt to its complex conjugate 
through 

5TL = ——-{2VLx5VL-TLTi[VtxVt5n]) (191) 



^The field of complex being kept distinct from the field of octonions, such that the Hermitian property 
is defined for the octonions only. 
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and one computes that Q extermizes Wa if and only if 

e-2^°det[n](2fi X DZ + Vi{Z -det[n]Ti [Vi'DZ])) + e^''^ (^ - detlilpZ) =0 (192) 

Because this equation is f/(l) x £'6(-78) invariant, one can solve it in any basis. There 
always exists an element g G -EgC-rs) and a phase e^^ G U{1) such that e^'^^Z and e'^^g{'DZ) 
are in the non-standard diagonal form defined in [26] 

e^'PZ = ex (e^° + ze"^" sin(2a))f? + e^-'"(ei + 6 + 6) 

/6 \ 



e^'^^(!DZ) = e-x (e-^° - ze*" sin(2a))^ 1 - e'^ 



06 
VO 6/ 



(193) 



In this basis one computes that (192) implies (up to a sign) 

e'^g{n) = e-'"-"^ 1 (194) 

and using then (188) that W = 2g. The 27 constraints (192) altogether with the positivity 
condition W > therefore determine Q in terms of Z and DZ, such that iy(r2) defines 
the single-centre non-BPS fake superpotential. 

Understanding the link between the potential Wa{^,Z) away from its extremum 
value Q^,{Z) is important in order to be able to generalise the triangular inequality for 
the central charge [5] to show that the energy of a composite non-BPS bound state is 
inferior to the sum of the energies of its non-BPS constituents. If the extrema of Wa (fi, Z) 
were necessarily maximums, it would follows directly that this is indeed the case since 

Wa{n, Zi + Z2) = Wain, z,) + Wa{n, z^) < Wa{n,{z^), z^) + 1^,(^.(^2), z^) (195) 

but the Hessian f^^ I „ /^-i can admit positive eigenvalues, and this identity is not trivially 
true. Nevertheless, there ought to be a similar property relying on the compatibility of 
the two constituents that would imply that this inequality holds. 

Extension to the interior of the solution 

We will now generalise these equations to describe the flow of the various flelds allover 
space. For this we deflne the spatial components of the vector fleld strengths as the 
horizontal components with respect to the connexion dt + u (with J = 1, 8 of S'f/(8)) 

Fij = dwij + Cijdu (196) 
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and define the functions of the scalar fields defined as the central charges acting on these 
vectors. Within our specific Ansatz (116), 

^^u 5MW] ~ 5MW] 

So. ^ 6t 

Using a complex basis for the 56 charges in terms of 28 complex Q^"^ and their complex 
conjugate Qij, one will have in general 

Cu = -e"''-^ - \e^"nuM (198) 

with 

nij = u'^jj]^ {u,, + n,,) - tH,ij\ {i^''' + ^'') (199) 

where we used the standard notation of [53] for the scalar fields such that 

% = Ui,''^Qu + Vi.uQ''^ (200) 

In the single centre case, TZu can be identified up to a renormalisation with the small 
vector appearing in [31] in the framework of A/" = 2 supergravity. Note that because 
the duality group acts linearly on the field strength, the direction in R^^ defined by the 
tensor TZu will always be a constant, as it is in the specific duality frame of the last 
section. Moreover, (180) implies that TZjj satisfies the constraints of a 1/2 BPS charge, 
and its stabiliser subgroup of E-j(j) is [54] 

^6(6) X E^^ C ^7(7) (201) 

The simplicity of the solution we described in the last section comes from the property 
that we parametrized the moduli in E'k^t) / SU^{8) in terms of the specific parabolic sub- 
group (R!J_ X -Eece)) ^ K-^^ which coincides with this stabiliser, such that the action of the 
scalars on TZjj amounted to a rescaling by the dilaton and 

- {ooij + Qij) = u,/-^nij + Vi.ijTZ'^ (202) 

was constant. In this way the defining tensors Uij and Qij are determined in terms of the 
scalar fields and the specific direction in R^® determined by the small charge TZjj of unit 
mass, noting that its norm is fixed by 

^z{n),,z{ny^ = i (203) 
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It is clear from (202) that in any other duahty frame, Uij and Qij will flow allover space 
accordingly. 

Using some identities valid when v'^v is a Hermitian matrix over the quaternions (no 
£ component)^ 

det[t;"ft;] : v^-^dJCv"^ : = Tr [{v^v) x {v^v)dJC] 1 - 2 : v{{v^v) x dK:)v^ : 

: v{HdY)v^ : = : v^~\i^{vv^) x {vv^)dY)v^^ (205) 

one computes 

e^Z{i.F) = l(dU- '-e^^ ^ du\ - ^Tr [v^-^dty-^] (206) 

and using the definition t = ti + it2, that 

e^'DZ{i<F) = -- (dU - -e"^ i. dco + -Tr [v^-^dtv]) 1 + -p(^) : t;t-M(tj - ft^)^-^ : 

(207) 
Therefore it follows that 

(it/ _ le^u ^^^^ If/ (^(^^^ _ rj.^ [®Z(*F)]) (208) 

which generalises (156). Similarly, the scalar fields momenta are defined as 

Dti - Wtz = -2ie^'DZ{^F) + ze^(Z(^F) + Tr ['DZ{i.F)]) 1 (209) 

Note that this equation shows that each split octonion component ^PQ+iyi^ of —2ie^'DZ{'kF) 
is equal to the '£ conjugate' of the corresponding split octonion component 'as — it/i' of 
Dti + iDt2. This implies that their respective contribution to the equations of motion of 
the scalar <;^ cancel precisely such that the latter are indeed constant. 

The property that Iie['DZ{-kF)] is defined as a 3 by 3 Hermitian matrix over the 
quaternions {i.e. that its i component vanishes) is also valid allover space. We conclude 
therefore that 

ooik^.i {Z^\F) + uj^m^^Zp^iF)) = (210) 



^To prove the second identity, one uses associativity to cancel the v factors on the left-hand-side, and 
one computes for the right-hand-side, using 6 = v\j for short 

6x6, , , 6x6 , , , , ,6(6x6) 6 x (ik{6 x 6)dY) 
-— r-r £^ /ix £Wy -— — = -2(6x Ui<(6x 6)dY)) ,J =-2 ^ ; ,^, ' '- 

1^6(^6 X 6)dY 



dct[6] 



li^dY (204) 
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also generalises allover space, such that in general one would have 

dU - -e^^ * rfw = --e^{VLijZ'^{i.F) - u'"' Zij{i.F)) (211) 

To compute the expression of the scalar fields momenta, we will use the property that the 
e(8) 0so*(16) momenta P is a nilpotent element characterised by a generator h of so*(16), 
which can be expressed in a Cartan basis as [.ofooooo ]• Defining the chiral Spin*{lQ) spinor 
P in terms of fermionic oscillators as in [37] 

\P) = (dU + -e^^ -kdu + e^ Zij{'kF)a'tt' H [ui/'^dvkuj - Vijijduu'^)a!'a^a!'a} 

+ -eijkipqrse^Z'^{^F)a!'a}aFa'^a'a' + -eijkipqrs{dU - ^e^^ ^rfwja^a^a'^a'aPaM |0) 

(212) 
one obtains that h can be chosen as 

h = — Vtija'a^ + -Vt'^aiaj (213) 

such that \P) = l-P'^') + |P'^') + |P'^') decomposes according to (167) and 

h \P) = |P(^') + 2|P('') + 3|P<'0 (214) 

Or, alternatively 

(h-l)(h-2)(h-3)|P) = (215) 

Solving this equation one gets back (210,211), and moreover that the scalar fields mo- 
menta satisfy 



Uij 



^"^dvkiij - Vijijduk/-^ = 6e^ (^uJiijQk\p^i]qZ^"i'*^F) - QiijZk^i-kF) 



- ^u[,,nki] (n^^zp'^i^F) - cu^'^Zp.i^F)) 

+^Ql,,Qki] {nP'^Zp.i^F) - Up,ZP''{^F))') (216) 

This expression is complex self-dual thanks to (210)/ The component in the (2®14j)]a of 
5^(1) X 5^(3) vanishes, exhibiting 'fiat directions' parametrizing the quaternionic space 



^One computes for instance that 



■^£ijklmnpq {'-^"^^^^^ ^'^'' Zrs — fl"^"^ Z^'^) — UJ[ijQk\p^l]qZ''"^ — ^[ijZkl] — - 

1 



+ -n^.j^ki] {nP'^Zpg - u^pgZ^") + 2f}[,,u;fcjj,f},],(Z^« + u;P™f7«"Z„„) (217) 
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F4(4)/(5'p(l) x^2 Sp{3)). Within our Ansatz Uij and fiij are constant allover the flow, 
as opposed to the phase a which appears in the BPS system as the triviahsation of the 
'modified Kahler connexion' [5, 55]. However, it is clear from equation (202) that this is 
an artefact of the specifically simple duality frame we chose to define the solution. 

This equation directly implies that there is always one trivial vector field in the 
system. Using the imaginary part of (210), one obtains 






(218) 



and therefore the vector field in the direction Darboux conjugate to TZjj is vertical, i.e. 
in C,{dt + u) only. 



4.2 Rotating attractor 

Let us consider the expression of the momentum P in the near horizon region. One 
computes that 



dU - -e^^ * du 
2 



XA 



1- 



iJ A cosf 



v/-/4-JicOs2(^^) 



) (dln(r) 



■J A sm{9A)d9j, 



2^-h-JlcosH 



(219) 



and similarly that 
(p(?):t;t-M(t,-a)t;-i:) 

Using these expressions one gets 
e^Z(*F)) 



1 - 



iJa cosf 



V-h - Jl 



cos^ 



J A sm{9A)d9A 
-I4 ~ J^ cos^( 



(220) 



2j^COs( 



2V ^-/4-j2cos2(e_,WV 



(d\n{ 



r] + I 



J A sm{9A)d9A 
^/-h-Jlcos\9A] 



e'^'DZii.F) 



1- 



2JaCOs( 



•^ 2 V ^-/4-j2cos2(e_,) 



-dHr)+^^d^^i£^]l 
^-h-Jlcos\9, 



(221) 



Let us introduce the following notation for simplicity 



w = (dU + -e^^ -k du 
V 2 



XA 



Z = (e^Z{i.F) 



XA 



S = J fp(0 : v^~'d{t, - et^)v~' ■) "DZ = (e'^'DZi^F) 

2 \ / XA 



XA 



(222) 
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We define the generator h (which can be expressed in a Cartan basis as [.iSooooo]) of 

50*(16) 



h w = e"2*^Z - Tr ['DZ] 
hZ = e^'^w - Tr [S] 



h 2)Z = e~^*"S - w 1 + S - ITr [S] 
hJ: = e^'^^Z-Zl + 'DZ-nT['DZ] (223) 

which according to [26] characterises the nilpotent orbit associated to non-BPS extremal 
black holes. It decomposes the coset component of e8(8) into representations of SU*{8) 
as 

128 ^ 1'-" © 28'"'^ © 70^'" © 28''^ © l'"' (224) 

Choosing the phase as 

gi^ ^ V-h - Jl cos^(^a) - JJa cos(g^) 

one computes that 

h Pr\xA = '^Pr\xA 5 h Pe\xA = ^Pe\xA (226) 

It appears therefore that Pr\xA li^s in the nilpotent orbit associated to a non-BPS extremal 
black hole, and Pe\xA li^s in the minimal nilpotent orbit. The outgoing momentum Pr\xA 
exhibits that the scalar fields take values determined by the local central charges. The 
pull backs of the field strengths on the horizon define the charges of the black hole, such 
that 

Fij\xA = QiJ sm{9A)d9A A d^A (227) 

It follows that the scalar fields (139) at 9a satisfy the attractor equation 



e^'^'DZ - Zt + 'DZ -ITi [l)Z] = (228) 

which reproduces the standard spherically symmetric non-BPS attractor equation, only 
at tJ7 = 0. In general it is associated to a symplectic matrix Qij of SU{8)/Sp{4:) such 
that 

^[ij^kl] + -^^ijklmnpq^"^"'^^'^ = (229) 

and the attractor equation reads 

^[ijZkl] + ■JTT^ijklmnpq^^^Z^'^ = (230) 

48 



The equation transforms covariantly with respect to the U{1) Ehlers symmetry, such that 
this extra phase is directly related to the phase of w. In our case 

a,j = e^'^Uij - Q,j (231) 

4.3 Non-BPS interactions 

To begin let us argue that the asymptotic tensors Uij and Qij associated to a solution can 
be determined in terms of the asymptotic moduli and the charges if the solution includes 
enough many black holes. In this subsection the moduli and associated central charges 
(as function of the scalars) will always be understood to be evaluated in the asjTiiptotic 
region. 

The set of charges Q^/ carried by the black holes will span a vector space R'^^ C R^^ 
which includes a Lagrangian subspace associated to the invariant symplectic form. We 
choose a Darboux basis associated to this Lagrangian, such that one can then define the 
orthogonal complement R^^ C R^^. In our duality frame, these charges are simply the 
ones admitting only as non-zero components p^ and a 3-vector of quaternions p. One 
can single out the p^ charge in an invariant way by using the property that the second 
derivative of the quartic invariant restricted to the adjoint representation evaluated for 
such charges is only zero if the 3-vector of quaternions p is null 



BQdQ 



133 



oc pp^ (232) 



The condition 



8^/4 



(7^) = is required for the vector TZij to be small [56]. The 

133 



dQdQ 

singled out direction determines therefore the vector 7?./j, and then the combination 



cu, 



«j 



+ Vtij through (202, 203). The obtained symplectic form Uij + Vtij decomposes in a 
single way into a sum of orthonormal rank 2 antisymmetric tensors as 

3 

ujij + ni^ = Y,^'ff (233) 

A=0 

Then using equation (210) on all the charges will determine u as 



a;,.- = uvr \ ^Qa Y1 '''^'^'ff {Z^\Qa) + uj^lio%,Z,,{QA)) = (234) 






This way one determines Uij and Vtij in terms of the charges Q^/ and the asymptotic 
moduli. We will therefore consider them as functions of the charges and the asymptotic 
moduli which transform as the central charges Zij with respect to Ej(j^ transformations. 
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It will turn out to be useful to define the real tensor in the 1 © 14 of 5*^(3) 

Wi^ = ^^ik^ji {z"' + n'^pn^'^Zpg) - ^n^jRe[u'''Zki] (235) 

This tensor is described within our Ansatz by the Jordan algebra element 
W{q, p) = Re[(DZ{q, p) - Z{q, p)l] 

= ^^^((^'[i^^n{^-^p')n-p{'^)--v'~\'P-ap')v~'-) 

= ^^ (Tr [{[ X [)p] 1 - pi,)v^-'pv-') (236) 

which is a Hermitian matrix over the quaternions because of (186). 

It is now time to have a second look at equation (135). First of all let us note that 
assuming that det[pyi] ^ 0, we can inverse this relation to 

1a[+ y{lA-lB)] : = TTTl + ^ ^22^) 

^A \xb-Xa\ det[pA] 

Then we will contract this equation with [ x [, and recognise each component as being 
writable in terms of the tensor function T4^ of the moduli. Indeed, one can define the 
following invariants by using the property that Tr ['Jt'"] is reproduced in general by the 
trace of Wik^^^ to the power n 
Tt[{[x[)pa] 



det[/'] 
Tr [[{pA X pb)] 



Tr [W{qA,PA)] (238) 

Tr ['M^(gA,PA)]Tr [W{qB,PB)] - 2Tt [W{qA,PA)'n^iqB,PB)] 



det[/'] 

'f I t1 ^ I 

8det[WiqA,PA)]+4^Tr [W{qA,PA)]Tr [WiqA^PA^] - ^^ ["n^iQA^PA)] 



det[pA] _o j„^r„,;^„ _ ^1 , ^rr„ [„,;/„ „ Mrr„ rr,.;^„ „ \2l rp^ [rrW^^ . ^ M3 



det[/'] 
For convenience we will rename these invariants as 

Tr[{[xC)pA] 



Wa = 

Wab 



det[r] 

Tr [[{pA X pb)] 

det[r] 



_ Tr [pAXpBPc] . . 

^^^^ = d^] ^^^^^ 

Using the quartic -^7(7) invariant Ii{q,p), we can also define the invariant quantity 

^ _ Tr['^(|t,-fj) 
Ja = 



^ .T. 



Waaa 

/I / Tfli 

{[ X C) [2{pA X Pa) X qA- ^^r [pa^A Pa + det[pA] [k. 



detM L^" -' V"^^" ^"^ ^"^ 2 ^^"^"^ ^" ' ^^^"^ V" det[r] 
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Finally, we can write the equations 



37^ 



m 



det[r] 



^{1a-1b)- 



W. 



B 



B=/=A 



Ju f 



Xa\ 



Ja_ 



(241) 



which permit to determine all 7^^ coefficients. It turns out that 7^ are not covariant 
quantities, but det[^]7^ — m are. One can define the A^ by A^ matrix 



M 



AB 



Cj^l \xc-xj_\ |a;2-a:i| 

Wi 3 + V ^C 



|x3-Xi| 

W3 



|xi-Z2| 



Cj^'i |a:c-a::2| |a;3-a;2| 



V 



(242) 



/ 



and solve 



det[/:]7^ -m = Y^ MaIJb (243) 

B 

We will now consider the symplectic products of (135) with pB, which give 

Tr [pa{pb X pb)] + Tr [psipA x pa)]' 



(7A-7ij)(Tr [C{pA X Pb)] 

^ v^ // xTr [pa{pb xpc)] 

+ 2^ ((7A-7c)- 



C^A,B 



X^ Xc 



ilB-lc) 



•^A ■^B I 

Tr [pa{.Pb X pc)] 



Xf 



Xr 



:Tr [^APs - ^bPa] 



(244) 



and can be rewritten in general 

/ \(^XT , ^abb + Wbaa\ 
{lA - 1b)[Wab + , , ] 



c^a,b 
where we used for short 



EA N WabC I N W^ABC 

I KlA - 1c)tz —\ " WB - 7c)- 

'-/- A n^ 



\X A Xr 



Xn Xr 



{Qa,Qb) 

det[C] 



(245) 



1 



(Qa, Qb) = ^{qoaPb - QobPa + Tr [(Ia'Pb - (Ib^a]) 



(246) 



Considering 7^ to be determined as (243), one has 26N remaining equations in (135) for 
a solution with A^ centres. Assuming all pA to be linearly independent, equations (245) 
provide inf[7V(A^ — 1)/2,26A^] covariant equations which constrain the centres in terms 
of the angles Uij, Qij. 

Note that, as opposed to the BPS case, one obtains inf [A^(iV - l)/2, 26A^] > 3iV - 6 + 
Sn,2 equations, and not only A^ — 1, such that the system does not even admit trivially a 
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solution for A^ > 10. In general these equations will give inf [A^(A^ — 1)/2, 2QN] polynomial 
equations of degree A^ — 1 for the A^ centres. This is because the leading terms of order 
N in the inverse radius all cancel. In the two centres case one gets 

'"^^ "^'i" (Ji-J2)i^i2-3(gi,g2) ^ ^ 

In the three centres case one can find the distance separations as real roots of quartic 
polynomials. Therefore they admit closed form formula in principle, although these 
expressions are not really illuminating and we will not display them. 

In general (for more than three centres which are not on the same axe, or more than 
four which are not on the same plan), equation (135) determines 15 A^ charges in terms 
of AN — 6 parameters, and therefore it is clear that the moduli space of such composite 
black hole solutions is only defined for a much smaller submanifold of allowed charges. In 
the axisymmetric situation, 15A^ charges are determined in terms of 2N — 1 parameters 
(respectively 3A^— 3 for the plan). It is therefore clear that considering the largest possible 
landscape of allowed charges requires to consider non-axisymmetric configurations, as 
opposed to the BPS case for which one could always rotate a configuration of charged 
black holes to an axisymmetric one. 

This property is confirmed by the appearance of a new contribution to the total 
angular momentum, when at least three centres are not aligned. In such a case, the 
charge factor reads 

Tfl 

7cTr [pA X pBVc] - -^^^^"^ ^PA X vbVc] 

-3Tr [^{qA,pA)\Ti [W{qB,PB)\^^ W{qc,Pc)] 
+4Tr ['M^(gA,PA)]Tr ['M/(gB,pij)'M^(gc,Pc)] 
+4IY ['M^(gB,PB)]Tr [W{qc,Pc)'n^{qA,PA)] 
+4Tr ['M^(gc,Pc)]Tr [W{qA,PA)'n^iqB,PB)] 

+8Tr [W{qA,PA)Ti^{qB,PBW{qc,Pc)]jJ2M~l,J{qn,PD) (248) 

^ D 

and because the second term does not contribute by symmetry, we obtain that the con- 
tribution to the angular momentum is a duality invariant function of the charges and the 
asymptotic moduli. 

In principle one would expect the expression of the angular momentum to not only 
be a duality invariant, but to moreover be independent of the asymptotic moduli such 
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that it could take an integer valued expression for arbitrary moduli. However, given the 
complexity of the complete expression of (165) in terms of the charges and the asymptotic 
moduli, and the property that the charges are themselves constrained, it is very hard to 
prove that this is indeed the case. The first step would be to actually understand clearly 
the constraints on the charges for the existence of a solution in terms of duality invariant 
equations. On the other hand, note that whenever one gets such contribution to the 
angular momentum, the different contributions to the angular momentum are then not 
aligned by definition, and there is no reason for the sum of the standard contributions to 
give rise to an integer valued total angular momentum anyway. 
Let us note nevertheless the following observation. We define 

dA = iAC+y^ [lA - 1b) I ^"^ I (249) 

such that (135) reads simply qA = 2pA x c{a- For generic charges pA, we know that the 
interactions disappear if and only if all 7^ are equal. Then one has ({a = tis since 

VA + VB \ , sr^ f, \ PC , \ Vc 



c^a-cCb = {'yA-lB)U + -, ^-^)+ X^ ((7a -7c) I ;-{lB-lc) 



(250) 
One checks straightforwardly that this implies that the charges are mutually local 

(Qa, Qb) = Tr [pA X pB{dA - de)] (251) 

which is the standard property for non-interacting black holes. Moreover, one can check 
that the sextic invariant proposed in [57, 58] to describe interacting solutions also vanishes 

dl dl \ 
^TT' ^7^) = 8det[pA] det[pB](det[i(4] - det[4]) ^ 2 det[pB\{qA{VA ^ {^a - d:B))qA) 

+ 2 det[pA] {q\{pB x (^a - 'iB))qB) - {qoB det[pA] + qoA det[;7B])Tr [pA x pB{dA - ds)] 
- 2[q\{pB X pB){dA - dB){pA X PA)qA) - 2{q\{pA x pa){(Ca - (Cb){pb x pB)qB) (252) 

when dA = ds- This suggests that this invariant might play a role in the description of 
these solutions. 



5 Almost BPS system 

Let us now discuss the generalisation of the almost BPS system of equations derived in 
[12] to A/" = 8 supergravity. 
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5.1 The solution Ansatz 

The system includes one real harmonic function V, a Jordan algebra valued (z.e. a 3 by 
3 Hermitian matrix over the quaternions) harmonic function /C, and a 3-vector over the 
quaternions of harmonic functions Y. It also includes a Jordan algebra valued sourced 
function £ satisfying 

di.{dC + JCxlCdV-Vd{ICxlC)) =0 (253) 

and a real function M sourced according to 

d^{dM -VTt CdlC+{Y^{dlC)Y)) =0 (254) 

It will be convenient to define 

= V det[C] - {Y^ C X CY) (255) 



and its derivatives 



and 



--^ = VCx C-2Cx (ryt) (256) 

^ VC~ {YY^) (257) 



dCxdC 
where iYY'^) is the Hermitian matrix over the quaternions built from the 3-vector Y . 

The metric is defined by the scaling factor 

e-4^ = - M^ (258) 

and the Kaluza-Klein vector satisfying 



The moduli are defined by 



92 

i.du = dM - Tr —r-^dlC (259) 

oC X dC 



K.zE±IL^^,e*'- (260) 



The Kaluza-Klein components of the vector fields are 

* dwo = Tr [CdIC - ICdC] + Vddet[JC] - det[JC]dV 
i^dzv = dC + ICx ICdV - Vd{IC x IC) + H {{dlC)Y - ICdY) 
i.dv=VdlC - ICdV - i^dY 
i,dv^ = -dV (261) 

Of course setting F = one gets back the system of A/" = 2 supergravity [12]. 
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5.2 Solving the differential equations 

Algebra shows directly that the poles of V and Y are incompatible with the poles of /C for 
the Noether charge carried by the pole to correspond to a regular orbit. The condition 
for regularity is the same and we will already assume that the function /C admits pole at 
the BPS centres only. This is in agreement with the result obtained in [13] for the STU 
model. For simplicity we will consider one single BPS centre at a; = 0, such that 

JC = l^+V{0)-^^ (262) 

whereas V and Y will admit poles at the (arbitrary Xa 7^ 0) non-BPS centres 

V = h-J2^^, Y = y~J2-^ (263) 

^ I A| A ' ^' 

This restriction will somehow be complementary of the case discussed in [14], where they 
considered one single non-BPS centre and arbitrary many BPS ones. 
One can then solve 

c = vior(h-y^^^Y-^ + ^^^±^^ + y-^ + [ (264) 

\ ^^^^^ ry \ ^ I nc ^ ry ^^^^^ nc ry 

A ' ' ' ' ' ' A ' ' 

The choice of parametrization is justified by the value of the electromagnetic charges at 
each centre. One computes the charges of the BPS black hole as 



g=^o = -Tr 



^ + 12a\^\\ 1 , , . , NT^ P\ 



^^*+(^x^ y^)v]-V{Q)-'det[p:,Y. 



(^ = q* + i*p, 



V{0) 
Y{0) 



I 1 3 



no) 

P° = (265) 

and the charges of the non-BPS black holes as 

qoA = -Tr [/C(x^)^^] + det[}C{xA)]PA 

0-^ = "^^ - IC{xa) X IC{xa)pa + ^^!^{xa)Pa 

'Pa = ^{xa)p\ + l-kpA 

Pa = Pa (266) 
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In order to compute the expression of the function M, we will need to define the 
function 

J An \x - y\\y\-^ \ dy'J \xA\\xB\\y - XaWv - XgWy]^ 

We will not compute this integral, however one can show that it converges for any value 
of X, such that the function Kab is regular everywhere, and satisfies the differential 
equation 

\^|3^a| \x Xb I |Xb||2; '•^b\) __ i 

Using this function one computes the solution for M 



AKab = V \ " ,, "' , " rr^ ■ V— (268) 



detfpl 1 /^-E4^^(l-F^) 

M=m+\/(0) — -— r h -ir [(^* + 2^ X ;7^)p^' 



+ 






X!(Tr [(/i^^ - PaOp*] + (2/Wa) + (PaP*2/))- 



2/ ,\ — v\ ± A J- j^ / r ^ J \€7 r '<'! ^^ / ' \i f\r "^ ij J I 

\JijA\\ •J-i \ \ JU JU A \ 

1 / X-X^~ 



I o 

/y» I ^ I ry I \ I rp ^ 

A A 



-^(O)^^E(^ATr [^^vA + (pWa)) , ^ PI , (l 
det[p*] v;^ , /3|x — x_4p + |x^P — 4|x^||x — x^l 1 1 \ 

^(0)"' Y. (P^Tr [-qj,pA +PBTr [?^F*] + (pWb) + {v^bV*Pa)) (F(A,b)o + , ,, ,, , 



+\/(0)"Met[p.] 5^ j)Vb fci? + 



X^ 3^ 



|2 



-B 



I "? I I "^ I I 

XO /y O rp 

A>B , A| I «| I I 

^0 



^ I SI ^ I -41 I I 

+ V^(O)- (/.Tr M - {y^P.y)) ^ - E 4rZT# (269) 



JU . JU JU A I 



where we chose the dipole harmonic function at x = such as to cancel the divergence 
in -^ in the scaling factor e~^^, and each pole such that u is globally defined. Fi^a,b)o 
is the function defined in the Appendix A for Xc = 0. We refer to [13, 14] for a more 
detailed discussion in the case of the STU model. 
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5.3 Near horizon geometry 

The BPS centre 

In the BPS centre near horizon geometry, one computes that 



k 



uj = ^^^^'^;';y^ + 0(1) (270) 

and 

_4U _ -q*odet[p^] + Tr [{q^ x q^){p^ x p^)] + {ql{p^ x p^)q^) - |(Tr [q*p*]y 

I 7 |2 .7 .^iv 2 

. ^ _det[^ >^ Pa 4 .979N 

an angular momentum induced by the interactions with the non-BPS black holes. The 
horizon area is thus 

S, = A7T^/h{q,p) + \J,\^ (273) 

The scalar fields turn out to do not sit at the attractor values, in fact they are defined 
at the horizon as (choosing coordinates such that Jldi = J^dz) 



with 



t = — -^ L. 

dq,o 

Indeed one can check that the BPS attractor equation is not satisfied 



(274) 



^Z ^ Jl + Jl sin^(g.) + I J* sin(g.) v//4 + P, + J,^ sin^(g.) ^ ^ 
2(J4 + J2 + J2sin2(^,))f 

although the scalar fields are entirely determined by the electromagnetic charges and the 
local angular momentum. Note that the solution is not quaternionic at x = (it admits 
non-trivial components in €) and therefore does not sit obviously in the A/" = 2 truncation 
in which Z{q,p) defines the central charge; all these components can nonetheless be 
eliminated by a duality transformation that shifts the axions by "^^y^m; such that it 
does. 

There is therefore no enhancement of supersymmetry in the near horizon geometry. 
The horizon is itself a squashed sphere with an angular momentum induced by the 
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interactions with the non-BPS back holes. As a consequence, the horizon area is not 
determined by the quartic invariant of its charges, but is increased by the norm squared 
of its induced local angular momentum. This horizon is larger than the BPS horizon in 
all directions, because its metric 

-Id -\~ \J* 



ds' = J U + \J,\^ + \J,\2sm\e,)de: + " ' ' =sin^(g,)V (276) 

verifies that 

T -L I 7' |2 

^^ + ''^*' > v^ (277) 

The existence of interactions seems therefore to create a squashed horizon that surounds 
the BPS horizon. 



Non-BPS centres 

The near horizon geometry of the non-BPS black holes is much simpler to extract. One 

directly obtains that 

J\{x — XaYcLx'' 

Jb "^ A I 

and 



u: = e.,, ^':_y.;-' + 0(1) (278) 



g-4C/ _ PAdet[gJ + (pJ4(g^XypA) _ / JAi{x-XAy \ 2 ^ ^ ^3 

Using the property that the charges of the non-BPS black holes are obtained by 'T- 
duality' of parameter /C(x^), one concludes that 

Sa = '^TTy/-h{qA,PA)-\JA\^ (280) 

as expected [40, 50, 51]. 

5.4 The asymptotic region 

The ADM mass 

To compute the ADM mass, it will be simpler to consider the asymptotic behaviour of 
the various functions in terms of the total electromagnetic charges 

V=h-^ + 0{r-^), Y = y-^ + 0{r"^) , /C = /t+ ^~^^ +0(r"^) (281) 
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and 

C= C+ ^ ^ — - — ^^ + 0(r-2) (282) 

The total go and q charges are determined such that 

go + Tr [^^]+Tr [(^^x^--^)p] - (^det[^] - ^Tr [^ ^])/ = 

q-k.p-Uv-k.p'')y = ^ (283) 

h 

Requiring the absence of total NUT charge implies 

M = m + ^-^ — -^^^-^ f^^^^-^ — ^ '^' + C r-2 284 

r 

This is indeed the case provided 



J- ^^^ / O \ T J- / \ *^ A "^ 



|2 



FiAM^) + I II I = Oir-') , Kab{x) + ,3, ;' = 0{r-') (285) 

in the asymptotic region. We prove explicitly the first limit in Appendix A, and we 
will assume that the second is valid. This is ensured by the conservation of the three- 
dimensional e8{8) current, since there is no poles outside the black hole horizons. 

By definition 

e-^^ = l + 4^^^ + 0(r-2) (286) 

r 

which implies 

h det[/'] = 1 + m^ + {y\[ x [)y) (287) 

Using this equation, the asymptotic scalar fields simplify to 

io = i+^^^(( X r - ^r X yyA +i^l (288) 

One computes that 



-'•^ADM , 



^ (- det[rK + /iTr [(r X 0(? + 2^ X p - ;t X ;t/)] 

- 2{y\ci + 2k.xp-k.y< I^P^)y) + {y\l x [)p + p\t x t)y) 

Off} \ 

+2mTr [/]tj9° - [p] + — (y^(p - ^P°)y)) (289) 
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After some algebra, one computes that it can be rewritten as 



M^oM = ^ ( Tr 



f X [ --[ xyy^)(l 



2Tr 



((t + <'*|)x(rxr-|fx,,t) 



^(rxf-|rxy;,t)^(fxr-|fx;,;,t,,y 



1 + m"^ 
2m 



det 
Tr 



l^^t--[xyy^ 



+ Tr 



k, + H^ X (^ + ^*|) X (/^ X r - ^r X yyt 



1 + m? 
where we used in particular that 



^ + ^*D "" {^^^-i^^yy^) X (^xr--rxyytjjjpo\ ^390) 



r X r - ^r X yy^) x(^[x[-^[x yy^^ = (det[r] - ^{y^[ x [)y)^ (r - ^) 



det 



1^ ^ ( - j^i: X yy^ 



det[[]-^{y^{fx[)y)"' 



(291) 



Before to be able to rewrite this expression in terms of the asymptotic central charges, 
we will first observe that 



N = : I go + Tr 



^+^*|)Q, 



Tr 



t+f*|)x(t+^*f) 



det 



/t + ^* 



^^^^^fx[-lcxyy^)x{[x[-l[xyy^]]'P 

y 



hi 



:Tr 



(x[-^[xyy^^ X (^[x[-^[xyy^^ x (^ + £*|)])/) (292) 



1 + m'^ 
vanishes identically according to (283). Then one can compute that 



-''^ADM — -f^^ADM 



3gia _^ g-3m 



-A^ 



^(3e*"Z(g,p) + e-3-Z(g,p)-e— (Tr ['DZ{q,p)] -Tr ['DZ{q,p)] 



4V 



-e^"Z(g,p) + --^- 



e -Im[e2*"Z(g,p) + Tr [®Z(g,p)]] 



for 



m + i 



yrr 



m^ 



(293) 
(294) 



This form indeed generalises the formula derived in [34] for the STU model. 
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Angular momentum 

Using the property that the angular momentum is defined by the asymptotic expression 
of the Kaluza-Klein vector 

00 = ^^i^-|'^^'^^' + o{x-') (295) 

one can compute that it also defines the asymptotic component of the function M 

M = m-J,^ + ^y^^ (Tr [{VC + h[)p,] - (Y^p^Y) - (y^p^y)) + Oix-') (296) 

where the last term includes the terms in A and -02 . To compute this expression, we 
would need the asymptotic behaviour of the function Kab- We will consider instead a 
undetermined parameter j\^ such that 

Kab = —, — r^, r^r^ - - ( 1 — T - 1 T ) 1 T " r^ ' JABiy^ + 0[x ') (297) 

The specific choice for the second term (that could have been absorbed in the constants 
Jab) ^i^^ become clear shortly. Using this expression one computes 



~-^' = wur^T^^u + 2i^^p*)p*]Yl 



2^(0)^^ Lv-z. '-■--.,,,, ,^|^^ 



7 

Pax\ 
2 



A 



^v{o)Y.{i^'' HhqA-PAOp*] + (^Wa) + ip\p*y))-^^ 

7n^J2(iPA^'^ ^'^aP*] + iPAP*PA)) 



detj^/ Y^ P% \ fY^ PAX\\ Y^ i 

^^TT^y Yl \ (Pa^"^ [?bP*] +PBTr [^^p*] + {PaP*Pb) + (pLp*P^)) 



yl>B 



-^ A -^ B I ^ ^ \ , -^ A "T -^S 



/y« rp \ rp /y I ly rp 

rf nn ■" ly '' ( rp np \ , ly ( ly " rp 

•^ a -^sI •'^ b y'^ a -^bJ •'^b\'^a 



rp rp rp I rp rp I rp I rp 1 

I A •'^B I I-'^aI |-^s| ^^l-^A •'^b\ ~r \-^A\ ~r \-^b\)' 

det[p*] >r^ // 1 1 \3x'4-x* 



2^(0)3 

^ ' A>B 
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After some algebra one obtains that the angular momentum can be rewritten as 

^ A 

+ 2 I] (^o^^B + Tr [a.A'P 

A>B 







Tr It , D ,-,1 n*^ n ,-, 1 ■* 




/ x^ - 


-Xb 


det[p*] ^-^ -j 





A>B 



A>B ^ 

I 19 7 / \ 7 I I 19 7 / \ 7 



rif rf* rvt 

. A -^B M^/l -if 



\\Xa ~ Xb\ + \Xa\ + ps|j 



(299) 



The first two lines correspond to the usual electromagnetic dipole angular momentum. 
Note that although the charges of the non-BPS black holes would have been trivially 
commuting in the absence of a BPS centre, the presence of the latter rotates the charges 
differently depending of their distance from the BPS centre. The third lines contains the 
angular momenta that can be felt in the near horizon region of each black hole, although 
note that the angular momentum induced by the non-BPS black holes in the near horizon 
region of the BPS one has the opposite sign with respect to the contribution to the angular 
momentum in the asymptotic region. The last line corresponds to a contribution to the 
angular momentum that would vanish if the centres were all aligned, similar as the one 
that exists for non-BPS black hole composites. Relying on duality invariance, it is clear 
that the constants j^^ should vanish for an axisymmetric solution, and that they should 
produce a contribution similar to the last line contribution in general. But we will not 
prove this property in this paper. 

6 Locally BPS system 

Let us now describe a locally BPS system, which is defined by the usual BPS system of 
differential equations with additional harmonic functions sourcing axions that do not lie 
in the M = 2 truncation. We consider therefore two scalar functions Lq and K^, and 
two Hermitian matrices over the quaternions C and /C, which altogether would define the 
harmonic functions of the maximal M = 2 truncation. In addition we will consider a 3- 
vector of quaternions of harmonic functions Y , such that only K^ , K, and Y are harmonic 
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dC 


- iJC 


dLo 


-iK^ 



functions, whereas Lq and C satisfy 

d ^ (dLo + Y\dlC)Y - Y^KdY - {dY^)lCY) = 

di<{dC + K^d{YY^) - YYUK^) = (300) 

The metric is then defined as in the BPS system by 
e-4^ = J4 = i^Odet[£] -Lodet[/C] + Tr [(£ x £)(/C x /C)] -^{K^U + Tr [JCC]^ (301) 

and 

i^du = - {K°dLo + Tr [ICdC] - Tr [CdIC] - LodK°) (302) 

which is well defined because 

K^ALo + Tr [JCAC] = K^Ti [ICA{YY^)] + Tr [IC{-K°A{YY^))] = (303) 

The scalar fields are defined as 

i^Y (304) 

and the three dimensional vector components of the vector fields are 

^ dwQ = Lo + Y^{dlC)Y - Y^ICdY - {dY^)ICY 
i.du> = dC + K°d{YY^) - YYUK^ + i^ {{diqY - KdY) 
i.dv = dlC + i^ {YdK^ - K°dY) 
*dv^ = dK^ (305) 

However, there is no regular composite black hole solutions which involve the functions 
y in a non-trivial manner, because a pole in Y necessarily renders the corresponding 
axion field singular. 

7 Maximal nilpotent orbits 

The three solvable systems of differential equations we discussed in this paper were all 
associated to nilpotent orbits which exist in D4, i.e. that they can be represented by 
elements of so(8, C) appropriately embedded inside eg. In fact the almost BPS system 
is associated to the principal orbit of 5*0(4,4), and the non-BPS composite (as well as 
the locally BPS) are associated to the subregular orbit of 5*0(4,4). They correspond 
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to two inequivalent nilpotent orbits in the coset component es 0so(16, C). Since they 
correspond to direct generahsations of the STU models, it was somehow natural that 
their expression was very similar to A/" = 2 supergravity expressions. 

In fact there are orbits of type -D5 which generalise these three orbits. The embedding 
of the principal nilpotent orbit of 50(5, 5) defines a system which generalises the almost 
BPS system, and the embedding of its subregular nilpotent orbit defines two inequivalent 
systems generalising the non-BPS composite system and the locally BPS system. Then 
the linearly realised SU{2) x SL{3, H) symmetry is broken to SU{2) x SU{2) x SL{2, H) 
and these three systems include 8 more functions which transform in the fundamen- 
tal of SL{2,M). The SO{5,5)/{SO{3,2) x 50(2,3)) model defines a truncation of the 
50(8, 5 + n)/(50(6, 2) x 50(2, 3 + n)) one that captures most information about A/" = 4 
supergravity black holes solutions. Indeed, the maximal nilpotent orbit of so (8, 5 + n) 
which admits a non-trivial intersection with the coset component is of type D^, and 
corresponds to the maximal nilpotent orbit of 50(5, 5). 

The maximal orbit of -^8(8) which admits an intersection with the coset component 
orthogonal to so*(16) is of type Eq, which justifies the use of the £'6{6)/5pc(8, R) to 
obtain black hole solutions in A/" = 8 supergravity. Once again the principal nilpotent 
orbit of Eq appropriately embedded inside e8(8) defines a generalisation of the almost 
BPS system, and the subregular orbit of Eq defines two inequivalent nilpotent orbits 
associated respectively to a generalisation of the composite non-BPS system and of the 
locally BPS system. In this case the linearly realised symmetry of these systems is broken 
to ni=o SU{2)i, and the latters include 4 more functions defining a quaternion, that is 
56 functions in whole. We will describe in this section the three nilpotent orbits of type 

Eq in C8(8)- 

It is useful to consider the decomposition of so* (16) 

3 

so*(16) ^ 0so*(4)i © 0(4, ® 4,) 

i=0 i>j 

3 3 

= sts" © su(2), © 2« ® 2^'^ © {2i © 2j) (306) 

with respect to which the coset component 128 decomposes as 

128 ^ 2^''^ © 2(^) © 2(^) © 2(^) © 2''> © 2^'^ © (2^ © 2^) © (2o © 2i © 22 © 23) (307) 

The four 5^2 correspond precisely to the STU truncation, as well as the 2'°* © 2*^^ © 
2(2) (g) 2(3) of the coset component. We will call Hi Cartan generators of the 5^2 ■ The 
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SU{2)i factors cannot be further decomposed within a real graded decomposition {i.e. 
with respect to a g[^ and not a u(l) which would require complexification). So for any 
graded decomposition, the (2o ® 2i ® 22 ® 23) will always be of grade zero, and the 
corresponding graded decomposition of C8(8) will accordingly always includes so(4, 4) in 
its grade zero component. 

This property implies that the nilpotent orbits of Cg in the coset component orthogonal 
to so(16, C) which admit a non-trivial intersection with the real 128 inside C8(8), are 
necessarily associated to weighted Dynkin diagram of so (16, C) for which the component 
associated to compact generators of the Cartan subalgebra of so* (16) are all null. They 
are of type [.„?„,„^, ] with c, 5, c, (f G IN [59]. 

We refer to [60, 61] for the complete classification of the nilpotent orbits of -£^8(8) • 

7.1 Almost BPS solvable algebra 

The maximal nilpotent orbit of -^8(8) which does admit a non-trivial intersection with 
the coset component is associated to the weighted Dynkin diagram of 63(8) [zoSozzz] which 
defines the graded decomposition 

28(8) = 1'""' © 1'-'°' © 1'"" © 8^-^' © (1 © 8a)'-" © (1 © 8e)'-'" © (1 © 1 © 8,)'-^^ 

© (1 © 8a © 8c)'"'' © (1 © 8„ © SaT^'' © (1 © 8a © 8^)'-'' © (1 © 1 © 8^ © 8e)'"'' 
© (gli © 0^1 © Sli © Qh © so(4, 4)) "" © (1 © 1 © 8^ © 8c)''' 

© (1 © 8a © 8,)'^' © (1 © 8, © 8a)'=*' © (1 © 8a © 8^)"' 

© (1 © 1 © 8c)"' © (1 © 8c)"' © (1 © 8a)'" © 8^"' © 1'^' © 1'"^' © 1'"' (308) 

A representative element of the orbit is a generic element of grade 1, for which all four 
components are non-zero and such that the vector and the spinor have a non-vanishing 
norm. The four components have all different weights with respect to the four q{^ such 
that no combination of them leaves the element invariant. The vector is only preserved by 
5'pm(3,4) C S'pm(4, 4), and the chiral spinor by 6*2(2) C S'pm(3,4). The corresponding 
nilpotent orbit has dimension 216, and is isomorphic to 

i?8(8)/(G2(2)X]R'^'+') (309) 

Its intersection with the coset component is associated to the so* (16) weighted Dynkin 
diagram [.^jo^o^o ] , {i.e. Hi + 2H2 + SH^ + 5Hq in the STU basis) which defines the graded 
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decomposition 

so*(16) ^ 1^-^°^ © (2o ® 23)^-'^ © (2o © 22)^-'^ © (1 © 2o © 2i)^-'^^ © (22 © 23)'"'^ 

ffi(lffi2i©23©2o©2i)'-*)©(2i©22©2o©22)<"''©(l©2i©23©2o©23)(-''©(2i©22©22©23)(-^^ 

3 

© (gli © sti © gli © 0^1 © su(2),) "" © (2i © 22 © 22 © 23)^^' 

© (1 © 2i © 23 © 2o © 23)'" © (2i © 22 © 2o © 22)'" © (1 © 2i © 23 © 2o © 2i)<'' 

© (22 © 23)'^' © (1 © 2o © 2i)(«' © (2o © 22)''' © (2o © 23)'«' © l'^"' (310) 

for which the 128 decomposes as 

128 ^ 1'-"' © l(-«' © (2i © 22)'-'' © (1 © 2i © 23)(-'' © (22 © 23)(-«' © (1 © 1 © 2o © 2i)(-^' 

©(2o®22©22©23)<-''>©(l©2o©23©2i©23)(-''©(2o©22©2i©22)'""©(l©l©2o©23©2o©2i)'-^' 

© (2o © 2i © 22 © 23)'°' © (1 © 1 © 2o © 23 © 2o © 2i)('' 

© (2o © 22 © 2i © 22)*" © (1 © 2o © 23 © 2i © 23)'''' © (2o © 22 © 22 © 23)('' 

© (1 © 1 © 2o © 2i)'^' © (22 © 23)(''' © (1 © 2i © 23)('' © (2i © 22)''' © 1<'' © 1<"' 

(311) 

where 2, © 2j is the real vector representation of S0{4:)ij = SU{2)i x^^ SU{2)j. Again 
an element of the orbit carries a non-trivial component in each irreducible representation 
of the grade one component of the 128. Therefore no combination of the q[i generators 
leaves it invariant, and only the diagonal of the three SU{2)i (for i = 0,1,3) leaves 
both the 2i © 22 and the 22 © 23 invariant. The corresponding orbit of Spin* (16) is of 
dimension 108, and is isomorphic to 

S pinl (16) /{SU (2) X SU{2) x R^+sj (312) 

Any element of the positive grade component is also an element of this orbit, as long as 
the grade one component admits generic elements, for every value of the higher order 
components. The components of positive grades define a solvable algebra of dimension 
108 which intersection with the coset component is of dimension 56, and therefore which 
permits to define a system of solvable equations for 56 sourced harmonic functions, i. e. 
one for each electromagnetic component. This system contains the almost BPS system. 

We will now describe this nilpotent algebra n in terms of the associated differential 
graded algebra /\n*, as we did in section 2.2. 
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Before to consider the maximal nilpotent orbit, let us restrict ourselves in a first 
step to the maximal nilpotent orbit of SO {4, 4) in e8(8) (discussed in more details in 
this paper), i.e. of weighted Dynkin diagram [oogoo^J [.4Sooo2ol- In that case the relevant 
solvable algebra inside e8(8) admits an £'6(6) automorphism and can be written in terms 
of the exceptional Jordan algebra. 

For the solvable algebra 

(1 © 27)^^' © 27^^' © 27'" © 1'*' © l^') (313) 

we define V, Eq and M associated to the grade 1, 4 and 5 singlets respectively, and K, E 
and L the grade 1, 2 and 3 Hermitian matrices. We chose the notations to be as similar 
as possible to the ones used in the paper such that one can recognise to which generators 
are associated the various functions of the associated solvable systems. The differential 
reads 

5E = \/K , 5L = E X K , 5^0 = Tr LK , 5M = EqV + Tr EL (314) 

It is nilpotent because the cross product x is symmetric. Decomposing these expressions 
in terms of split real and split octonion numbers according to 



/ K^ k, k* 



K 



3 J^2 



k* K^ ki I (315) 

V k2 k* K'' 

and analogously for E = {E\ei) and L = (Lj, P), one gets 



6E'-- 


= VK' , 6ei = V\^i ■ 




6Lr- 


= E'+^K'+^ + E'+^K'+^-e 


jk^ + kjG^ , 


sr-- 


= -E^k, + jre, + e*^_2k*+i- 


k* p* 


6Eo-- 


= >;(L,ir^ + rk*-ka"), 





SM = EoV + J2{E'Li + eir-re*) . (316) 

i 

where i = 1,2,3 and repeated indices are not summed. The nilpotency can be checked 
explicitly by using the property that the real part of a square is symmetric, and that 
the real part of a cube is cyclic. To decompose the algebra in its so* (16) and coset 
components, one must decompose the split octonions in terms of quaternions 

J^i ^ rCj + tSi , Gj ^ ej + It/i , Ij ^ tj + iXi [oil) 
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and use the Cayley product (40). 

The algebra then decomposes in terms of quaternions as 

5U = K+^K^^^ + K+^K'+^ - e,kl + he* + y,s* - s.y* , 

5x' = -E'si + K^Vi + k*^iyi+2 + ki+2yi+i - ej+2Si+i - e*^^Si+2 ; 
SEq = J2{LiK' + rk* - hr* - x's* + Six") ; 

i 

SM = EoV + Y,{E'Li + e,r* - Te* + x'y* - Vix'*) . (318) 

i 

which disentangles the elements V, K^ , ki, yi, Li, P , M associated to the coset component 
from the elements Si,E\ei,Xi,Eo of the so* (16) subalgebra. The automorphism group 
then reduces to SU{2) x 5L(3, H). 

One obtains the maximal algebra by adding three generators associated to split oc- 
tonions Tj = Tj + iti. ti will be generators in the coset component 128 whereas rj will 
correspond to generators of so* (16). The E^^q) automorphism will then be reduced to 
Spin{4:, 4), such that octonions only get multiplied through the triality invariant products 
which preserve Spin{4, 4) in the three trial representations pi ^ 

(p,((xy)*))* = p,+i(x)p,+2(y) (319) 

According to this rule, the grading and the Jacobi identity determines completely the 
algebra up to a redefinition of r^, as 

6E' = VK' , 6ks = K^s , 6r2 = -r^r^ ; 
6K' = kgr* - rgk* , Sk^ = Kh^ - r^k* , ^63 = Vk, + Eh^ ; 
6E' = VK^ - rge* + egr* , 6k, = Kh, + k*r; - r*k; , Se^ = Vk^ + Eh^ - r*e* ; 
6K'^ = kir* - rik* + ksr^ - rsk^ , 5L3 = E^K^ + E^K^ - egk* + k^e^ , 
6ei = Vki + Eh, + e^r; - r*e; ; 
5E'^ = VK^ - rie* + eir* - rge^ + e2r; , 51^ = -E^k, + A'^ei + e^k; - k*e; - Lgri ; 



®Here pi, p2, ps state for the same element of Spin{A, 4) in the chiral spinor, the antichiral spinor and 
the vector representation, respectively. 
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SL2 = E^'K' + E'K^ - esk; + k^e^ - lir* + nl* , 

61^ = -E^k2 + K^e2 + elk; - k*e; - Lgra - l^r^ ; 

61" = -E^kg + K^es + elkl - k^e*, - L2I3 + r^l^* - l^^r* ; 
6Li=E^K'^ + E^K^ - eik* + kie* - Fr^ + rsl^* - Fr* + rgl^* ; 
6Eo = Y,{L^K' + rk*-k,r*) ; 

i 

6M = EoV + J2{E'Li + eiV* -Ye*) . (320) 

i 

Note that there is a particular truncation of this algebra for which there is an enhanced 
50(5,5) symmetry, namely for r3 = 0. The algebra then corresponds to the nilpotent 
orbit [00S0022J) [AooioI- a representative element of the maximal nilpotent orbit in the 
coset component has V,K^,ts,ti all different from zero. 

7.2 Locally BPS solvable algebra 

There is another pertinent orbit which also defines a system involving 56 harmonic func- 
tions, and which contains the BPS system. It can be obtained from the nilpotent orbits 
of e8(8) weighted Dynkin diagram [^ogozozl which defines the graded decomposition 



= 1'-^' © 2*-'' © 8'-'' © (2 © 8a)<-'' © (1 © 2 ® 8cy-'^ © (2 © 8a © 8J'-'' 
© (1 © 2 ® 8a © 8^)<-'' © (2 © 2 ® 8^ © 8c)'-'' © (gl^ © Qi^ © gl^ © sk © 5o(4, 4)) '"' 
© (2 © 2 © 8„ © 8,)'^' © (1 © 2 © 8a © 8,,)'" © (2 © 8^ © 8^^^ 

© (1 © 2 © 8e)<'' © (2 © 8„)''' © 8f © 2*'' © l'**' (321) 

Such a nilpotent element admits a non-trivial component in each irreducible representa- 
tion of the grade one component, including an element of the 2, a doublet of two linearly 
independent non null vectors and a non null chiral spinor. The two orbits associated to 
this graded decompositions are distinguished by the property that the two vectors are 
of the same signature or not. If they are, the doublet of vectors is left invariant by an 
SU{2, 2) C Spin{4:, 4), and the chiral spinor by SU{2, 1) C SU{2, 2). Otherwise, they are 
left invariant by an >S'L(4) C S'pm(4, 4), and the chiral spinor by 5*17(3) C SL{4). These 
two orbits are of dimension 214, and are isomorphic to 

E8(8)/(5f/(2, 1) K E^^) , Esis)/{SL{3) x R^'^) (322) 

respectively. There is one intersection of each of these orbits with the coset component. 
The intersection of the second is associated to the weighted Dynkin diagram 



■ 2202020 
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Its positive grade components do not include the 8 charges of the STU truncation, and 
we will only discuss it in the following subsection. The intersection of the former is 
associated to the weighted Dynkin diagram [.0502040 1 (^-e. Hi + 2H2 + 4if3) which defines 
the graded decomposition 

3 



;i©2©2o©22©22©23)(-''©(2(g)2o®2i©2i(8)22)*""©(g[i © gli © s^i ©s^a © 0su(2)i)' 

i=0 

© (2 © 2o ® 2i © 2i ® 22)''' © (1 © 2 © 2o © 22 © 22 © 23)*'^ © (2i © 22 © 2i © 2^^<'> 



(0) 

i=0 

1 © 2 © 2o © 23)'^* © (2i © 23)<^) © (22 © 23)'«> © 1^'^ (323) 



with respect to which the coset component decomposes as 

128^2(-''©(2o©2i)(-«)©(2©2o©22)^-''©(2©2i©22)'~''©(2©2o©23ffi2o©22)^-'^ 
© (2 © 2i © 23 © 2o © 2i)(-'' © (2 © 2 © 22 © 23 © 2o © 23)*-'^ © (2o © 2i © 22 © 23)™ 
© (2 © 2 © 22 © 23 © 2o © 23)''' © (2 © 2i © 23 © 2o © 2i)<'' © (2 © 2o © 23 © 2o © 22)*'' 

© (2 © 2i © 22)*'' © (2 © 2o © 22)'" © (2o © 2i)('" © 2''' (324) 

A representative of the nilpotent orbit is then defined as an element of the grade one com- 
ponent with a non-zero 2 element, a doublet of linearly independent vectors of 50(4)23 
and a vector of 6*0(4)03. The latter is only left invariant by the diagonal 5'f/(2)o3 
subgroup of 50(4)03, and the doublet of vectors is only left invariant by a subgroup 
50(2) X 50(2) C 5L(2) x 5t/(2)2 x 5f/(2)o3. The element in the 2 is not left invariant 
by the 50(2) factor, such that the orbit is isomorphic to 

5pm*(16)/(f/(2) kR9) (325) 

and is of dimension 107. The associated system involves 56 harmonic functions and 
includes the BPS system. The property that one needs a doublet of linearly independent 
quaternions shows that such element does not exist inside the real e6(6), and therefore 
the -E'6(6)/'S'Pc(8, K) truncation might not be enough to describe this system of equations. 

Similarly as for the almost BPS system, let us consider in a first place the -£^6(6) 
invariant -D4 type nilpotent algebra associated to the orbit [oogoo^o] [.2S00020 ]• In that case 
one simply replaces the generator V associated to a D6 by a V associated to a D6, which 
gives 

5K = VE , 5L = E X K , (5M = Tr EL , 6E0 = MV + Tr LK (326) 
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or after decomposition 

5K' = VE' , 5ki = Ve, 6si = Vyi 
5Li = E'+^K'+^ + e^+2k^+i _ e^k* + he* + y,s* - s,y* 

Sr = e*^2k*i+l - K+2e*i+l - E'^i + K'^i - Si+Wi+2 + yi+is*+2 
Sx' = K'yi + k*^iyi+2 + ki+2yi+i - E\si - ej+sSj+i - e-^^Sj+a 
6M = Y,{E'Li + e«^" - ^'< + x'y* - y^x'*) 

i 

SEo = MV + J^i^i^^^ + ^^^i ~ ^*^" ~ ^'^*i + ^^^'*) (22'') 

i 

The maximal solvable algebra can be obtained by adding again the same three octonions 
Fj. And because the corresponding generators commute with the generator associated to 
V (as can be checked from the grading), one gets the algebra 

5^2 = _r3e* + egr* , Sk^ = Vk^ + Kh^ , 5r2 = -r*r; , ^62 = -r*e* ; 
6K^ = VE^ + kar; - rgk* , 6L3 = E^K^ - egk^ + kge^ , 
5ei = Ehi + e*r; - r^e^ , ^ka = Ve2 + K^Y2 - r^k^ ; 
5E^ = -rie* + eir* - rge^ + esr^ , 61^ = K^e^ + e^k* - k^e^ - L^vi , 
6ki = Vei + Khi + k^r; - r^k; ; 
6K^ = VE'' + kirl - nkl + kar; - rak^ , 
6L2 = E^K' + E'K^ - e^kl + kse^ - Iir^ + ril* , 

61^ = -E^k2 + K^e2 + elk; - k^e^ - Lgrs - l^r^ ; 

61" = -E^ks + ir^eg + e^kt - k^e^ - Lsrg + r^l^* - I'^r^ ; 
6L1 = E'^K^ + E'^K^ - eik* + kie* - ih; + rsl^* - ih; + rgF* , 
6M = E^L2 + E'L^ + ^ (e,r * - Fe*) ; 

i 

6Eo = MV + XI (^»^' + I'k* - ka") . (328) 

i 

Where we have ordered the terms in function of the grading, and we set E^ to zero because 
it is of grade zero. Note that a closer look to the algebra permits to recognise the SL{2) 
automorphism of the maximal BPS system (generated by E^ and its conjugate). 
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7.3 Non BPS solvable algebra 

Now let us come back to the orbit associated to the weighted Dynkin diagram 
{i.e. —^Hq + '^Hi + |if2 + IH3), for which so*(16) decomposes as 

50*(16) ^ 1'-^' © (22 ® 23)(-«' © (1 © 2i ® 23)(-^' © (2i © 22 © 2o © 23)'-*^ 

©(l©2o©23©2o©22)<"''©(2o©2i©2i©23©2o©22)'"''©(l©2o©2i©2i©22©22©23)(-'' 

3 
© (0^1 © 0^1 © gli © 0^1 © su(2)i) '"^ © (1 © 2o © 2i © 2i © 22 © 22 (8) 23)^'^ 

j=0 

© (2o © 2i © 2i © 23 © 2o © 22)''^ © (1 © 2o © 23 © 2o © 22)^'' 

© (2i © 22 © 2o © 23)*"' © (1 © 2i © 23)''' © (22 © 23)*'' © 1*'' (329) 

and respectively does the 128 

128 ^ 1<-'' © 1'-'^ © (2o © 2i)'-«' © (1 © 2o © 22)'-'' © (1 © 2o © 23 © 2i © 22)'"*' 
©(l©2i©22©2i©23)<""©(l©22©23©2o©22©2i©23)<-"©(l©22©23©2o©23©2o©2i)'-'' 
© (1 © 2o © 2i © 22 © 23 © 1)"" © (1 © 22 © 23 © 2o © 23 © 2o © 2i)<'' 

© (1 © 22 © 23 © 2o © 22 © 2i © 23)''' © (1 © 2i © 22 © 2i © 23)'=*^ 
© (1 © 2o © 23 © 2i © 22)<'' © (1 © 2o © 22)''' © (2o © 2i)(«' © I*'' © l^^^ (330) 

An element of the orbit is defined as a generic element of grade one. The four components 
are not left invariant by any combination of the qI^^ and the three 2j © 2j are only left 
invariant by the diagonal SU{2) of the four SU{2)i, such that the orbit is of dimension 
107 and isomorphic to 

Spin:{16)/{SUi2) K E=^^=^+^) (331) 

The positive grade component defines a solvable algebra which almost contains the one 
associated to non-BPS solutions in the STU model, but a singlet generator associated to 
one D2 is in the grade zero component. However, adding such an element to a general 
element of grade one does not modify the stabiliser subgroup. Let us write L*"', L^"' the 
coefficients of the elements (l©2j©2j)'"' in so* (16), and X,Xij the elements of the grade 
one component of the 128. We will consider the element of the 2j © 2j as 2 x 2 matrix 
which multiply themselves through the contraction of their indices associated to the same 
SU{2). Note that any 5*0(4) vector is invertible as a 2 x 2 matrix. One computes that 
the elements of the positive grade component of so*(16) which leave invariant a generic 
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element of the grade one component of the 128 are defined from 10 parameters as 

-^23 ■ ^23 = , L{^ = Xo3 XqiL , Lq3 ■ X03 = , L^2 = -^03 -^23-^03 ; 

r(2) Y _n r(2) _ v-'^y r(2) r(2) _ v-ly -l y y r (2) Ci-iO^ 

-t^Ol'^Ol — U, -L^3 — A Ao3-t^oi ' -^02"^ ^01 ^23-^03-t>oi (^OOZj 

An element which includes moreover a non-trivial element of grade associated to a D2 
with coefficient Y is also trivially invariant with respect to the diagonal SU{2), and is 
also left invariant by the positive grade elements as long as 

It follows that one can associate the nilpotent element to the deformed grading associated 
to 

~Ho + ^i/i + ^H, + \h,+ ^- [H, + H, + H, + Hs) (334) 

which gives the following graded decomposition of 50* (16) 

so*(16) ^ l*-'-^' © (22 ® 23)(-^-^' © (1 © 2i ® 23)'-^-^' © (2i © 22)'-*-^' © (2o © 23)^"^^ 

©(l©2o©23)'-=*-^'©(2o©22)<"''©(2o©22)'"'"^''©(2o©2i©2i©23)<-''©(2o©2i)(-'-^) 

3 
© (2i © 22 © 22 © 23)<''' © 1'-'+'-' © (gli © gli © Qi^ © qI, © 0su(2)i)'"' © l^^-^^ 

j=0 

© (2i © 22 © 22 © 23)(^> © (2o © 2i)f^+^) © (2o © 2i © 2i © 23)'"^ © (2o © 22)^'+^' © (2o © 22)^'^ 
©(l©2o©23)'^+^'©(2o©23)(^'©(2i©22)<"+^'©(l©2i©23)<^+^'©(22©23)'«+^'©l(^+=) 

(335) 

and of the 128 

128 ^ 1^-*-^' © l'-'-'^' © (2o © 2i)(-'^-^' © (2o © 22)'^'-^^ © 1<-'' © (1 © 2o © 23)^-^-^) 

© (2i © 22)*-^' © (2i © 22)*"'-^' © (1 © 2i © 23)'''^ © (1 © 2i © 23)^-'-^^ 

© (22 © 23 © 2o © 22)'"" © (22 © 23)'-'-^' © (1 © 2o © 23 © 2o © 2i)<"'> © l'"^' 

© (2o © 2i © 22 © 23)"" © 1(^' © (1 © 2o © 23 © 2o © 2i)('' © (22 © 23)''+^^ 

© (22 © 23 © 2o © 22)''' © (1 © 2i © 23)'^+=' © (1 © 2i © 23)'^' © (2i © 22)''+^' 

© (2i © 22)*"' © (1 © 2o © 23)*^+^' © 1'^' © (2o © 22)*'+^' © (2o © 2i)(«+^' © 1^'+^^' © 1'^+^' 

(336) 

Although this is slightly more subtle in this case, one gets in this way a solvable algebra 
which defines a solvable system of differential equations involving 56 harmonic functions 
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which generahses the composite non-BPS system of the STU model. Of course we could 
fixe e to any fixed value, for instance e = 1 permits to exhibit an additional SL(2) 
symmetry of the system. The graded decomposition is then defined with respect to 
2Hi + 3H2 + 47^3, and gives 



50*(16) = l'-**' © (22 ® 23)'-^' © (1 © 2i ® 23)'-''' © (2i © 22)<~'^ 

,(0) 



(1 © 2 © 2o © 23)'-*^ © (2 © 2o © 22)'"'' © (2 © 2o © 2i © 2i © 23^'"'' 



©(2i© 22 ©22 ©23)'-'' © (gli©0[i©0[i©sl2©0su(2)i)' 

j=0 

© (2i © 22 © 22 © 23)('' © (2 © 2o © 2i © 2i © 23)'" © (2 © 2o © 22)^" 
© (1 © 2 © 2o © 23)'*' © (2i © 22)''' © (1 © 2i © 23)'^' © (22 © 23)(^' © I'**' (337) 

and of the 128 

128 ^ 2^-^' © (2o © 2i)f-^' © (2o © 22)^"'^ © (2 © 2o © 23)'"'' © (2 © 2i © 22)'"*' 
© (2 © 2 © 2i © 2sY~''^ © (2 © 22 © 23 © 2o © 22)'"" © (2 © 2o © 23 © 2o © 2i)^-'' 
©(2o © 2i © 22 © 23)"''©(2©2o©23©2o©2i)'''©(2©22©23©2o©22)'"©(2©2©2i©23)('' 
© (2 © 2i © 22)'" © (2 © 2o © 23)('^' © (2o © 22)''' © (2o © 2i)('' © 2*^' (338) 

In order to obtain the 'maximal' non-BPS nilpotent algebra from the almost EPS 
one, one can again simply replace the K generators associated to D4's by conjugates 
K associated to D4. Using the already derived commutation relations, it remains only 
to compute the ones involving r^ and K. But noting that the algebra admits an SL(2) 
automorphism with respect to which (K, L) transforms as a doublet and r^ is invariant, 
one gets these relations without effort. 

5r2 = -r^r^ , ^k^ = -KsVi , 6\' = -LsV^ ; 
6K2 = -kir* + rik* , 5L2 = -hrl + nl^ , 

6k^ = -K^V2 - K^l , 6\^ = -Lsr2 - l^r* ; 
6E' = K^L-i + i^3L2 - k^i* + l^k^* , 

(5k3 = -A'2r3 + r;ki*-k2*r* , (51^ = -L2r3 + r^l^* - l2*r* ; 
5Ki = -khl + r2k2* - k^r* + r3k^* , 6L1 = -\hl + r2l^* - I^r* + r3l3* , 

^63 = -K^l^ + L-sk^ + k^n^* - F^k^* + Ehs ; 
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SE"" = K^Li + KiL-i - k^P* + l^k^* - rse* + egr* , 

5e2 = -K2\^ + Lak^ + k^*!^* - \^*]<i^* + E^r2 - r*e; ; 

5ei = -Kil^ + Lik^ + k^*!^* - l=^*k2* + Ehi + e*r; - r*e; ; 
6E^ = K1L2 + K2L1 - kH^* + Fk^* - ne* + eir*, - rae^ + e^r^ ; 

6V = J2{E'Ki + eik'* - k'e*) , 6M = ^(E% + e,r* - Fe*) . (339) 

i i 

The role of the three new quaternionic functions Tj in the system is clearly to relax the 
constraint (210), and to allow the scalar fields t to carry a non-trivial I component in its 
imaginary part. The Ansatz is significantly more complicated in this case. Moreover the 
system of differential equations is also much more involved, and the functions V and M 
are now sourced by polynomials of order seven in the harmonic functions. The definition 
and the analysis of the regularity of these solutions requires further studies. 
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A Solution of the Laplace equation with three point 
sources 

We define the solution to the Laplace equation 

^Fa,bc = I 2(^-^^)-(^-^^) (340) 

as the integral 

P [ dh {y-XB)iy-Xc) , . 

^A,BC = - / -^-: n M 131 13 '^'^^^> 

J 2tv \y- x\\y - XA\\y - XB[^\y - xd^ 
This integral is convergent for any value of x. Being regular everywhere, this integral is 
determined to be the solution of the Laplace equation which admits no poles and which 
vanishes in the asymptotic region x — ?■ 00. 
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As such, it follows that if Xa,Xb and Xc are aligned, the function is given by the 
axisymmetric solution 

pAxsym _ ^ / I -^ -^A | \-^ B -i- A\ PC -^^A | 

A nc^ — / \/ \li II II II II II I 

(342) 
We note also that by definition the function is invariant with respect to the exchange of 
X and Xai therefore whenever x is aligned with respect to Xb and Xc, it reduces to 

1/ rr> rr> rf ry If 'T' 

/ •^*^A *^S'^ r^C*^ 

i'A.BC 



( rr> rf 1,1 'T' 'T' I \ If 'T* 'T* 'T' O^ 'T' If If If 'f If If 

yuj Jjb) \^'^ -^c) \ l-^A •'^b\\-^a •'^C] \-^b •'^cW-'^a ■'^b\ I'^'b •''cII'^a -^c 

(343) 
The symmetry with respect to the interchange x -H- x^ will play an important role in the 
following. 

However, if the axisymmetric solution is not too hard to compute, the general solution 
is rather difficult to obtain. We will study the asymptotic expansion of the general 
solution in the asymptotic region in this appendix. But let us first show that the integral 
(343) indeed converges for all values of x. To see this, let us study the limit ?/ — )■ Xg in 
the most dangerous case, i.e. when x — > Xg as well. 

Limit X — )■ Xs 

We consider the integral on a very small ball Eg of radius e surrounding x^ , such that x 
lies in the ball. In this case we can expand the integrand in y — x^, (note that |x — x^l < e) 

_ f A {y-XB)iy-xc) ^ ^,. 

Jbb 27r \y - x\\y - XaWv - XbI^Iv - Xc\^ 

(X - X„)-(Xn - X^) f\ f^ , n ^OS 6* 

ar a cos u 



nr' rp \ \ nf nf "J T* 


Xq 1 


yX X]^)'yXQ ^b) 




rf rf -T* 'y "J -T* 


Xb\ 


[X — XsJ'yXc — Xq) 





J-1 \/r'^ — 2cos6'r|x — x^l + |x — x^P 

l^-^^l , 2r r , 2|x-Xs| 

dr— + / ar- 



"JF ^b\ J\x-xb\ "^'^ 

I II 131 iV^-^^f^) ^344) 

rf rf rf rf iJ rf rf \ \ -Ir^ I 

Let us now compute the explicit expression of the function Fa,bc at x — ;■ x^. In order to 
do this computation we will first decompose R^ into the interior of the ball Bb surrounding 
X in the limit |x — x^l << e, and the exterior of this ball. In the interior one gets the 
same result as in (344), with — )■ 0. In the exterior, one can expand the integrand 
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in the harmonics centred a.t x = Xg, such that 

d^y (y - XB)-{y - Xc) 



r^\Bb 27r \y - x\\y - XaWv - XbI^Iv - Xc\^ 



I -pL fa-^'Hy-^-) +o(,_,,) (345) 

Jn\Bn 27r \y-XA\\y-XB\'^\y-Xc\-^ 



d^y {y-XB)-{y-Xc) 

'R\Bb 

The singular component in e vanishes by symmetry, and one can therefore define the 
hmit at X — )■ Xs of Fa,bc as 

1- p {x-XB)ixc-XB) f d^y {y-XB)iy-Xc) . . 

lim l'A,BC = ] n 131 1 ~ / "o~ I n 141 13 ^'^^^^ 

x^xB \xa- Xb\\xc - Xb\'^\x - Xb\ Jji3 271 \y -x^Wy -XB^ly ~Xc\-^ 

Up to a harmonic function in x^, one obtains 

[X Xb)'[Xc Xb) 
^A,BC 






+ 1 ^ 121 II ^-1 n r^ + 0{x-XB) (347) 

rf i rf rf rf rf ^ rf rf rf rf ^ 

I A ■'^b\ \-^ a -^cW-^b ■'^c\ \-'^ a •^cW-'^b -^cI 



Using the property that this expression agrees with the axisymmetric one (342), one 
concludes that the possible additional harmonic function in Xa must vanish when Xa is 
on the line {xbiXc)- Moreover the integral is regular in the limit Xa — > C)0, and one 
concludes that (347) is indeed the right expression. 

Note that using the symmetry with respect to the interchange x f-)- x^, one obtains 
as well the limit of Fa,bc as Xa — > Xg 

J-, \Xa Xb)'{Xq Xb) 

^A,BC 



rf rf rf rf rf rf O 



+ ^^^l^^^^^^r^^2 - ] i^ ^ + 0(x.-x,) (348) 

rf Zr rf rf rf rf ^ rf rf rf rf Zr 



Limit X ^ X, 



The limit x — ?■ x^ is clearly given by the convergent integral 

IP . ^ I d^y {y-XB)iy-xc) 

^A,Bc(Xa) 



27r \y - XA^ly - XBl-^ly - Xc\ 
d d f d^y 1 



Oxb dxc J 2n \y - XA\'^\y - XB\\y - Xc\ 
d 9 \ f d^y 1 



dxA dxc ' """^J J 2tt \y - XA\^\y - XB\\y - Xc\ 
d^y {y - XA)-iy - Xc) 2 



ztt \y XBiiy Xa\ \y Xc| \Xa Xq\ \Xb X(^ 
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(349) 



and therefore 

A^^Fa,bc{xa) = -4] 7T, n7 (350) 



We conclude that 



Fa,bc(x.) = -2-. ^^^ ■'^y^^c Xa) ^ ^^^ 

for a harmonic function in x^, symmetric in the interchange of Xg and Xq- Comparing 
with the axisymmetric case (342), one obtains that Hbc must vanish whenever Xa.Xb, Xc 
are ahgned. Moreover, analysing the asymptotic behaviour of Fa,bc{xa) in the limit 
\xb — Xa\ — !■ oo, one obtains that Hbc = 0{\xb — Xa\~^)- We conclude therefore that 
Hbc must vanish identically and that 

Fa,bcM = -2i (^--■'^y(^o-XA) ^352) 

Limit X — )■ oo 

Let us now study the limit a; — )> oo. We consider a ball Sa centred at of radius 
A >> sup(|xa|, I^bI, |xc|), but such that \x\ » A. We decompose the integral into the 
integrals over the interior and the exterior of the ball. On the exterior, y is very large 
everywhere and one can expand the integrand as a Laurent series in y 

d^y {y-Xsyiy-Xc) 



r3\Ba 27r \y - x\\y - XaWv - XB\'^\y - Xc\'^ 
d^y [ 1 {xA + '2xB + 2xc)-y 



R3\B^ 27r \\y~x\\yY> \y-x\\y 

1 , {Xa + ^Xb + 2Xc)-X , ^^^^^_5x 



(^(r^; 



+ ^^ + ^(^ 



3|xP 5|xp 

1/1 (xa + 2xs + 2xc)-a; 



0(x"=') + C(A^3) (353) 



A^ \|x| 3|x|^ 

On the interior, x is very large compared to y itself, and one can expend the integral in 
the harmonics centred at x = 0. 



d^y {y-XB)iy 



Jb C 



Ba 27r \y ~ x\\y - XA\\y - XB\^\y ~ Xc\^ 
1 f d^y {y - XB)-iy - Xc) 



\x\Jb. 2n \y - XA\\y - XB\^\y - Xc\^ 



X f d^y^ {y-XB)iy-Xc 



13 / 9 y-l II 131 |3 +^(^ ) (354) 

x|^ Jb. 271 \y~XA\\y-XB\'^\y-Xc\-^ 



78 



By definition the terms in A cancel each others, as one can straightforwardly check at 
this order, by computing the contribution to these integrals in the neighbourhood of the 
sphere \y\ = A. Defining La.Ux-^oo{F) as the Laurent series in x of the expansion of the 
function F at a; — )■ oo, we conclude that 

Lau^^oo(^^,Bc) = - lim / -^:—-. rLaUy^oo i n Rl R 

A^oo Jj^3\Ba 27r li/ - xl \\y-XA\\y-XB\^\y-Xc\-^ J 

- I ^Lau._ (^\ , '''7'^-'^'"-' , (355) 

Jr3 27r \\y-x\J \y - XA\\y - XBV\y - XcV 



Note that the first term is also a Laurent series in x, because each term in the Laurent 
expansion in y will have a definite scaling in y, such that the corresponding integral will 
have a definite scaling in x (after having taken the limit A — )■ oo) . 
Applying the Laplace operator, one directly obtains that 

ALau,^oo(F^,5c) = lim / d^y5^^\x - y)Lau,^o. ( . %-^^)-(^-^c) \ 

^^^Jr^xba \\y - XAily - xbI-^Iv - xc\y 

= Lau_(^fcf^iH£_l£oL\ (35e) 

\ rfi rtr> rp /y O ly ly "-* / 

\iX/ iX/y^J/ Jb q\ \Jb Jb Q \ J 

because 

Lau,_,oo(5^'Ka;-2/)) =0 (357) 

and Bk does not contain x (before to take the limit) by definition. Therefore it follows 
that LaUa;->oo [Pa.bc) satisfies the correct Laplace equation as a Laurent series. 

In order to confirm that this formula makes sense, let us also consider the asymptotic 
behaviour of F^^bc in the limit x^ — )■ oo. F^bc is by definition symmetric with respect 
to the interchange x -H- Xa, and therefore the same formula must hold in this limit, i.e. 

T (17 \ V f ^^y ^ J ( {y-^BHy-xc) \ 

LaUa;^^oo ( -f^4,BC ) = - lim / —-- rLau^^oo i n isi n 

A^oo Jj^3\B^ 27r ly-x^l \\y-x\\y-XBV\y-XcV J 

d^yj / 1 A {y - XB)-iy - xc) 



Lau^4_^oo -, 7 -, n ni n (^^S) 

E3 271 - Wy-x^lJ \y - x\\y - XBl'^ly - Xcl'^ 

Applying the Laplace operator, we obtain the same result, but now only the second term 
contributes 

3„A(3)^. _ „U .„ f^_\ ^iy-^B)iy-Xa) 



ALemx^^^(FA,Bc) = / d^y6^''\x - y)LaU:i.^^oo ( i \ ] 

Jrs \\y-XA\j 



ly-XsPly-XcP 



T / _^_\ 2{x-XB)ix-Xc) , . 

Laua^-A^oo 1 f n m iT (359) 
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It follows that the two terms are necessary for Lau^^oo {Fa,bc) to satisfy the correct 
Laplace equation and to be the asymptotic series of a symmetric function in x and x^- 

Using the asymptotic expansion (353,354), one therefore concludes that 

d^y {y - XB)iy - Xc) 



Fa 



BC 



\x\ J 2-n \y - XA\\y - XBr\y - XcY" 

-^./j^„ (y--My-^c) ,o(,-3) (360, 

|x|-^ J 2t\ Xy-XAWy-XBYXy-XcY 

where the integrals are now computed over R^. These integrals are both regular at 
any value of x^, and it follows that they define the associated solutions of the Laplace 
equations in x^ which have no pole in x^, as we just discussed in general for the Laurent 
series. For instance 



d^y {y-XB)iy-Xc) _ 2(x - Xb)-(x - Xc ) 
2% \y — x\\y — XB\'^\y — Xc\'^ \x — Xb\'^\x — Xc\ 






and therefore 

^ y [y •^Bj'Ky xc) \Xb xc\ \xa Xb\ \Xj^ x^ 



2'Tr n I rf \ \ 1 1 IP \ '^ \ 1 1 T" " If T* If If If T* 

" \y -^AWiJ -^B] \y -^Cl \-^B •^CM-^A -^bW-^A -^C] 

And similarly 

d^y {y -XB)-{y- Xc) _ 2{x - Xb)-{x - X c) 

2-n \y — x\\y — XB\^\y — Xcl"^ \x — Xb\^\x — Xc\ 



(362) 



-^ I — y ... -.IL. ^ |3U,_^ 13 = ^ U_^ |3U_^ 13 (363) 



implies that 



d^y {y-XB)iy-xc) 
27r \y-XA\\y-XB\^\y-Xc\^ 

X 4 + Xn + X(-. i / Xq Xc 



rr> rr> rr> ryt ryt \ /y /y rf rr> 

I A •^SlI'^A -^d l-^S -^Cl \ l-'^A -^d l-^A -^B] 

c) ps ^ Xc\ (X^ — X^j — [Xb — Xc)'\Xa — Xb)\^b ~ Xej ('ip,A\ 

rf rp ly ly /y /y I /y /y _\_ /y rr> _\_ rr> rr> I 

\-'^ B ■^CW-'^A •'^bW-^A •^C\\\-^A -^sI ~r \-^A •'^C\ ' \-'^ B •'^ C \ ) 

In order to find this solution one first observes that 

-X4 + Xg + ^c Z\X Xb)'\X Xc) 



A 



XII I I I "^ I I "^ 
/y /y /y> /y /y O ly n/i O 



\Xb Xc\ (X Xg) [Xb Xc)'\X Xsjipi-B X^j /nnr\ 

+ 2 , 7T, r5 (305) 

rp ly O ly nr* O 
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The second equation 



AG 



BC- 



I I "^ I I "^ 



= 2[\Xb -XcI^V^s - (Xb -Xe)(xB -Xc)-V^j,g)- r pr (366) 

\ / HP 'T* 'T* 'T* "J 

can be solved by considering in a first step the equation 

AFbc = 1 1 R = J ^ (367) 

Expanding the harmonic factor, one computes that a particular solution to this equation 
is obtained as the integral 

da 

'BC ■ 



a/p^ + {z - ZcY^ p^ + {z- Zc- a{zB - Zc)y 



In' ^"^" 



jVZEfzfil! 



z-zc+^/p^+{z-zc)^ 



{Zc - Zb)\/p^ + {z- Zcf 

1 ( {x-xb)-{xo-xb)+\x-xb\\xc-xb\ 
y (x-a;c)-(a;c-a;B)+|a;-a::cl|2:c-a;B| 

nn nn nn 



(368) 



The corresponding Xg derivative only involves the derivative with respect to x — x^ 
because of the projection and gives 



^/ l^B ^c\ \^ ^b) [_^B ^QJ-yX Xb)\^b ^c) 

^Bc = \ n n \ 

1 



rfi I rf rf I rf rf rf rf I rf rf rf rf 



(369) 



However the second factor is singular whenever x goes to the segment [xb , Xc] , whereas 
Gbc is regular everywhere. We consider therefore the harmonic function 

YT Ps ^ Xc\ (X — Xg] — [Xb — Xc)'[X — XBJyy^B ~ Xcj /■o^n\ 

^Bc = , n M 1 13 /^Uj 

rf rf rf rf rf rf 

\'^ B ^O \ \^ ^ B \ \^ ^ C \ 

1 1 

X 



rf I rf rf I rf rf rf rf I rf rf rf rf 

\zb - Za\pe^y{z - ZbY + P^ + ^J{z - Zgf + p^) 

^{Z - ZbY + pV(^ - ^C? + p^{^ {z-zbY+P^^ {z-zcY+p^+{z-zb){z-zc)+p'') 



which interpolates between the singular harmonic function Vlog(p^) = ^^— near the 
segment [xs,Xc], and — — — 3- in the asymptotic region. Adding this contribution, one 
obtains the everywhere regular function 

^Bc-^-, n M ill m m rr ^'^'^) 

B ^ C \\ ^ B\\^ ^ C \ \ \ "^ B\ ~I~ "^ ^ C \ ~^ \'^ B ^ C \ I 

We have therefore the expansion 

I B -^d H^A -^Sl \-^A -^d 



F, 



A,BC 



W^B '^CW^A ^B\\^A **^C 

X / -X^ + Xs + Xc 1 



XO \ I n/-> rp I I rp rf I I rp /y I I I rp rp I I rp rp 

I \|iX/4 dU fj \ \dU y\ Jb (^ \ \JUj:^ Jb (^ \ \ \Ji/ j^ Jb (^ \ \Jb /i^ Jb 



+ 0{\x\-^) (372) 



This result is consistent with the axisymmetric solution, because the third term vanishes 
identically whenever Xa (respectively x) is aligned with Xb and Xc- One can also check 
that in the limit Xa — ?■ x^, this expression further reduces to 



^XjX Xg)-[XQ Xg) i \ Z' X ■ Xf 

XI II 1*^ I I9/I1III 
rf nc nf O nf ry» \ i^ I \ ry rf 

[Xa Xb)'\Xc -^b) J- 



'A,BC=\-, n fT ^ 1 [^ I I 77 -^ I ,f ) +0{\x\ ^,Xa-Xb) 

\ I 'T' 'T' I I 'T' 'T' I '^ I 'T' O^ I ■" / \ I 'T* I I 'T' I "^ / 



rp rp rp rp rp rp O rp rp rp rp ^ 

Jj J/^iX/y^ Jj^Jj^ Jb Q \ \Jb Jb (^ \\Jb Q Jb (^ \ 

^X_^Oa^XB) ^ ,^ ^^ _ ^ 

rp \0\rp rp \0 

Jy \JbQ Jb Q \ 

which coincides with (348) in the limit x — )■ 00. 

This confirms the asymptotic expansion of Fa,bc- For the composite black hole 
solution discussed in the paper, we are in fact interested in the function symmetrised in 
A,B 

F{A,B)c + 1 rr^ n 1 (374) 

^ ' rp rp rp rp rp rp 

\-'^A ■^'dl-'^fl -^CW-^ •'^C] 

rp I rp rp rp rp /If O" O^ 

Jb I Jj j\ Jj Q Jb Q Jb j\ ^ Jb (~; Jb J^ Jb Q 



2bP 



rp rp rp rp rp rp rp rp rp rp rp 

A Jj -D \ \ Jj t:/ Uj r^ \ \Uj T3 ijb A \ IUj A Jb f'^ \ \^ A "^^**-'R -Xj r 



A *^B||*^B '^C\ \'^B ^AW'^A ^C\ \^A *^CM**^S '^ C \ 

rp " I rp rp \ rp ( rp rp \ , ( rp rp \ I rp rp \ rp 

I A ^B\ \^B -^C) ^ \-^A -^B) \^B -^^Cjy-^A ^B) -^ 

i"o~i i~i i~i T/"! j j j j Tx" ~i 

rp O rp rp rp rp rp rp { rp rp I rp rp I rp rp I 

X\ \Xa Xb\\Xa Xc\\Xb Xc\y\XA Xb\ ~r \Xa Xc\ T^ \Xb Xc\) 
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B Bases for the solvable subalgebras 

In order to deal with nilpotent orbits representatives, it will is convenient to decompose 
the general charge ^ G 42 = e6(6) 0sp(8, R) according to a Cartan basis of sp(8, R). 
Up to redefinitions of the harmonic functions that define the solutions, one can restrict 
oneself to a Cartan basis in the purely non-compact component 20 according to the 
decomposition sp(8, R) = u(4) © 20, i.e. for which only Xabc is non-zero in (14). The 
most general such quartet of generators is parametrized by f/(4), and we will restrict 
ourselves to a particular case for simplicity. There is Cartan basis which is naturally 
associated to the D0-D4 configuration in the STU model, and which is defined as the 
generators Hq, Hi for which Xabc is defined as in (35) with only Qq and P* non-zero, and 
respectively equal to 



Ho:Qo = P' 



1 

'2 ' 



Hi:Qo = P' 



1 pj 

2' -* 



h for J ¥" ^ 



(375) 



The elements of the coset component decompose with respect to their weight in terms of 
these four generators. For the charges which sit in the STU truncation, one has 



DO : 


;i, 1,1,1) 


A^ = 


-s. = -g. = -p' 






D2 : 


;i, -1,1,1) 


M = 


— Si = S2 = S3 = 2-P 








;i, 1,-1,1) 


M = 


Si = — S2 = S3 = 2-P 








;i, 1,1,-1) 


M = 


Si = S2 = -S3 = |p2 






DA : 


;i, 1,-1,-1) 


N = 


—Si = S2 = S3 = gi = — g2 = 


-Q3 


= P0 




:i, -1,1,-1) 


N = 


'^'i = —^2 = ^3 = —Ql = Q2 = 


-Q3 


= po 




;i, -1,-1,1) 


N = 


'^'i = '='2 = —^3 = —Ql = —Q2 


= Q3 


= p° 


DQ : 


;i, -1,-1,-1) 


M = 


— Sj = ^Qo 






DO : 


-1,-1,-1,-1) 


N = 


-E, = Q, = P° 






D2 : 


-1,1,-1,-1) 


M = 


— Si = S2 = S3 = — 2-P 








;-i, -1,1,-1) 


M = 


Si = — S2 = S3 = — 2-P 








-1,-1,-1,1) 


M = 


Si = S2 = -S3 = -ip^ 






DA : 


-1,-1,1,1) 


N = 


—Si = ^2 = 23 = — gi = g2 = 


Q3 = 


-po 




-1,1,-1,1) 


N = 


Si = — H2 = S3 = gi = — g2 = 


Q3 = 


_pO 




-1,1,1,-1) 


N = 


'^'i = '='2 = —'='3 = gi = g2 = - 


-Q3 = 


_pO 


DQ : 


-1,1,1,1) 


M = 


— Sj = —oQo 







(376) 
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where all the unspecified charges are zero. Out of the the STU truncation one has 



D2 : 


(1,1,0,0) 


-v' 


= v =y = -pi 




(1,0,1,0) 


Ci 


= -(3 = Z= -P2 




(1,0,0,1) 


-e 


= e = K.=-Pz 


DA : 


(1,-1,0,0) 


Vo 


= vi = y = -<i^ 




(1,0,-1,0) 


c 


= e = -^ = -q' 




(1,0,0,-1) 


eo 


= ^3 = t = -q' 


Y : 


(0,0,1,1) 


-^0 


= rii = y = q^ 




(0,1,0,1) 


-c° 


= e = -^ = q' 




(0,1,1,0) 


-^0 


= ^3 = K. = q' 


T : 


(0,0,-1,1) 


v' 


= r]^ = -y = -pi 




(0,-1,0,1) 


Ci 


= (3 = Z=P2 




(0,-1,1,0) 


e 


= e = -K. = -P3 


fiat : 


(0,0,0,0) 


•Ji 




D2 : 


(-1,-1,0,0) 


-v' 


= v^ = -y = Pi 




(-1,0,-1,0) 


Ci 


= -Cs = -Z = P2 




(-1,0,0,-1) 


-e 


= e = -^=P3 


DA : 


(-1,1,0,0) 


Vo 


= vi = -y = Q^ 




(-1,0,1,0) 


c 


= e = -z = q' 




(-1,0,0,1) 


Co 


= i3 = --^ = q' 


F : 


(0,0,-1,-1) 


Vo 


= -vi = y = q^ 




(0,-1,0,-1) 


c° 


= -e = z = q' 




(0,-1,-1,0) 


Co 


= -^3 = -^ = q' 


T : 


(0,0,1,-1) 


v' 


= v^ = y = pi 




(0,1,0,-1) 


Ci 


= Cs = -^ = -P2 




(0,1,-1,0) 


e 


= e = ^=P3 



(377) 



Each system of differential equations associated to a given nilpotent orbit is defined 
by possibly sourced harmonic functions valued in the positive eigenvector of a Cartan 
generator. For the BPS system, the Cartan generator is simply Hq, and the system is 
generated by all DO, D2, DA, D6 charges, which all carry eigenvalue 1. The single-centred 
non-BPS system is associated to h{Ho+J2i Hi), and is generated by the charges D6, Y, D2 
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and -DO, carrying respectively eigenvalue 1 and 2. The almost BPS system is associated 
to 2Ho + J2i Hi, and is generated by the charges D6, DA of eigenvalue 1, F of eigenvalue 
2, D2 of eigenvalue 3 and DO of eigenvalue 5. The non-BPS system is associated to 
YliiHi, and is generated by the charges DA,D2 of eigenvalue 1, Y of eigenvalue 2, and 
-D6, -DO of eigenvalue 3. 

Another interesting system is the one associated to the maximal nilpotent orbit, for 
which the Cartan generator is 5-/^o + -^i + 2H2 + 3-^3. It extends the almost BPS system 
by bringing the charges T into the game, giving rise to a system that includes as many 
harmonic functions as there are electromagnetic charges 

20 = D'0 + 6xD)2 + 6xD)4 + D6 + 3xF + 3xT (378) 

Similarly, the maximal locally BPS system is defined by the Cartan generator H2 + 2H3 + 
AHq and includes 

20 = D)0 + 6xD)2 + 6xD)4 + D'6 + 3xF + 3xT (379) 

and the maximal composite non-BPS system is defined by the Cartan generator 2Hi + 
3H2 + AH3 and includes 

2O = D)0 + 6xD2 + 6xlM + I56 + 3xy + 3xT (380) 
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